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Abstract 

We study, using Rindler coordinates, the quantization of a charged scalar field inter- 
acting with a constant (Poincare invariant), external, electric field in (1+1) dimensionnal 
flatspace: our main motivation is pedagogy. We illustrate in this framework the equiva- 
lence between various approaches to field quantization commonly used in the framework 
of curved backgrounds. First we establish the expression of the Schwinger vacuum decay 
rate, using the operator formalism. Then we rederive it in the framework of the Feynman 
path integral method. Our analysis reinforces the conjecture which identifies the zero 
winding sector of the Minkowski propagator with the Rindler propagator. Moreover we 
compute the expression of the Unruh's modes that allow to make connection between 
Minkowskian and Rindlerian quantization scheme by purely algebraic relations. We use 
these modes to study the physics of a charged two level detector moving in an electric 
field whose transitions are due to the exchange of charged quanta. In the limit where 
the Schwinger pair production mechanism of the exchanged quanta becomes negligible we 
recover the Boltzman equilibrium ratio for the population of the levels of the detector. 
Finally we explicitly show how the detector can be taken as the large mass and charge 
limit of an interacting fields system. 
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1 Introduction 



Shortly after Hawking's discovery of black hole evaporation jl], Unruh || showed that a uni- 
formly accelerated detector moving in flat space perceives the Minkowski vacuum to be ther- 
mally populated at temperature TJy = a/2n. On the other hand, Heisenberg and Euler || 
showed in 1936 that the vacuum state of a charged quantum field interacting with a static 
electric field is unstable and decays into pairs. In 1951 Schwinger j|]gave the expression of this 
decay rate, using the technique known today as the Schwinger proper time representation of 
the functional integral. One can use the same formalism to describe these two phenomena, 
wherein the second occurs in the secure framework of usual quantum field in flat space. In this 
way we have a physically relevant model whose interpretation is unambiguous and which allows 
to exemplify several formal developments and check their validity. In this paper we mainly 
study, using Rindler coordinates on the (1+1) dimensional Minkowski space, the quantization 
of a charged field interacting with a background, constant, electric field E.The main motivation 
of this work is pedagogical: we hope to illustrate (and check) in the framework of an exactly 
(but non trivial) solvable model several conjectures that often are taken for granted in quantum 
field theory on curved spaces. In the framework of a charged quantum field interacting with an 
external constant electric field, we explicitly show that, the formal evaluation of vacuum decay 
rate based on an expression of (0, out\0, in) as a functional integral on quantum fields leads to 
results in accord with standard calculations,even in Rindler coordinates, where the integrals are 
no more gaussian. This supports the use of a similar approach for more complicated problems, 
for which no exact solution is known. We also illustrate several aspects of physics that we 
believe to be important; for instance: 1) the occurrence of a boundary effect, whose analog in 
the framework of the physics of the Boulware vacuum in Schwarzschild geometry erases the 
Hawking radiation; 2) the validity of the thermodynamical equilibrium relation in (and only 
in) the limit of heavy systems; 3) the significance of the "rindlerian energy" balance 4) the 
contribution of the Schwinger-like and Unruh-like vacuum fluctuations on the transitions of the 
detector, similar to the Schwinger and Hawking mechanisms of charged black-holes evaporation. 
Our paper is organized as follows. Section 1 is just a summary of the main body of the text. In 
section 2 we recall the definition of Rindler coordinates and classify the classical trajectories of 
charged particles in a constant electric field. The results of this analysis are intensively exploited 
when we interpret the subsequent results. In section 3 we recall the derivation of the Schwinger 
effect in the framework of the usual quantization (using the creation and annihilation operators 
formalism in Lorentzian coordinates) of a charged scalar field interacting with an external con- 
stant electric field. Then, in section 4, we perform the same work in Rindler coordinates. The 
results differ from those obtained using global coordinates by boundary terms. We interpret 
this difference as the manifestation of the same mechanism leading to the difference between 
the Boulware's and Unruh's vacua in the framework of black hole physics. Then, in section 5, 
we reconsider the problem in the light of the functional path integral method, and illustrate the 
validity of the standard approximation scheme (saddle point - W.K.B. approximation) used to 
evaluate functional integrals, by comparing their predictions to the results obtained previously; 
we show that both methods (fortunately) agree up to boundary terms. In the same way, we also 
reinforce the interpretation of the usual minkowskian propagator as a sum over winding Rind- 
lerian propagators || . Indeed using an expression for what we believed to be the exact Rindler 
propagator to compute the rate of particle creation from the vacuum, we shall recover at one 
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and the same time its leading term and the boundary correction terms. We also compute, in 
section 6, the Rindlerian particle content of the Minkowskian vacua. In section 7, we study the 
behavior of accelerated charged detectors. We show that the detector transition (in the limit 
where the charge of the exchanged quanta vanishes) are dominated by the induced transition 
of the detector with preexisting vacuum fluctuations that maximally overlap the detector world 
line. We also establish precisely how to interpret a two level charged detector as a limiting 
case of large mass M and charge Q, but of finite acceleration Q E/M, of an interacting field 
H |7| system. Finally, some mathematical appendices provide the technicalities underlying the 
results of the main text. 



2 Classical trajectories in Rindler coordinates 

We recall that Rindler coordinates divide Minkowski space into four patches, denoted hereafter 
by the labels R(ight), L(eft), P(ast) and F(uture). On each of these patches the coordinate 
transformations between Minkowski (t, z) and Rindler (r, £) coordinates are given by 



On each of these quadrants, the respective Rindler coordinates run from — oo to oo. In the 
following, from time to time, we shall also make use of another standard Rindler coordinate : 



which is obviously always positive. 

Let us emphasize that on R and L the vector field d T is timelike, pointing respectively to the 
future and the past while on F and P it is d% that provides us timelike directions, pointing 
toward the future on F and the past on P. In the following, we will be led to consider the 
various components of the boundaries of these patches. In two dimensions, the boundary of the 
full Minkowski space consists into four points: the future and past timelike infinities, denoted 
by i + and i~, the spacelike left and right infinities, i° L and i° R , and the four components of the 
null infinity: X R , T£, X R and X£. Moreover, the light cone issuing from the origin O defines the 
acceleration horizons of stationary Rindler's observers, whose trajectory equations are £ = C te 
on R, L and r = C te on P and F . It splits into four branches (constituted by the null rays 
emerging from its vertex O or ending on it): H^, H.\ and Ti^. They constitute parts of 
the boundary of the different Rindler's patches. The other components of their "boundaries" 
Qare the points i + , i~ , i° L , i° R and the eight pieces of the null infinity components denoted 

1 Between quotation marks because they are not really boundaries of regions in Minlowski space but bound- 
aries of the conformal Carter-Penrose diagram. In the rest of the paper, we shall not make explicit this 
distinction. 






(2.1) 
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by T^p (the future null infinity boundary of the R quadrant), 1p R , Ilfi ■ ■ ■ > ^-pr ■ All 
these geometrical considerations are summarized on fig. (1) which represents the well-known 
conformal Carter-Penrose diagram of (1+1) dimensional Minkowski space-time. 

A constant electric field, which in flat coordinates is given by JF = Edz A dt with E = Cst, 
reads as Q T = tEe 2a ^d^ A dr with e = +1 on R, L and e = — 1 on P and F. This field is 
invariant with respect to the Poincare group, acting on the (1+1) Minkowski space. It can be 
derived from the potential: 



E E 
A = —(zdt — tdz) = e ~^~ e 



2ai a~ l dr 



(2.2) 



Classically, it accelerates positively charged particles towards the right (z > 0), negatively 
charged particles towards the left (z < 0). The Hamiltonian describing these classical motions 
is 
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(2.3) 



the momentum p T is a constant of motion: 
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(2.4) 



and, once we impose the evolution parameter to be the proper timef] along the trajectory, the 
mass shell value of the Hamiltonian is fixed : H c i = — m/2, which determines p^. 
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(2.5) 



On R and L, where e = 1, f is of constant sign on timelike trajectories (f 2 = £ 2 + e 2a ^), while 
p^(= +m£e 2a ^) vanishes at the turning points £ + and £_: 



1± (1 - 2ft) 



1/2' 



(2.6) 



when Q = uoqE/m 2 a is less than ~. At these turning points: 



qE 
2amQ 



1 + (1 - 20) 



1/2' 



(2.7) 



showing that, for Q < 0, f_ is negative on the trajectories for which ^ < £_.The interpretation 
of these latter is obtained by noticing that particles moving forward in time (with respect to 
the global time orientation of Minkowski space), move with f > on R and f < on L, while 
anti-particles go in the opposite way. 
Finally, the acceleration at the turning points is: 



'( = ±ae~ 2a ^(l-2tt) 



1/2 



(2.1 



2 When no confusion arises we do not make explicit on which patch the Rindler's coordinates are considered. 
3 A dot indicating a derivative with respect to it. 



4 



Classical trajectories are branches of hyperbola (see fig. (2)). The geometrical significance of 
the constant Q is obtained by noticing that : 

is the squared invariant distance from the origin O to the center of the hyperbolic trajectory 
(see fig. (2)). 

On R and L, according the values of Q, we distinguish three classes of trajectories. 

• Those, on which > |. They have no turning points and correspond to motion with £ 
varying between from — oo to +00, with £ > 0. The centers of these hyperbolas lie in the 
sector P or F. Moreover, f > 0, so these trajectories describe particles in the R sector 
and anti-particles in L; more precisely, they describe: 

— On R, particles [1] (the numbers between brackets refer to the trajectories of fig. (1)) 
entering from P (with £ > 0), across the past horizon 7-Q with £ > and running 
towards the null infinity component or [2] coming from (with £ < 0) and 
leaving R via Tl^, 

— On L, anti-particles (particles going backward in time) moving between the infinity 
and horizon components I^f an d with £ < 0, or between Z^ P and7i^ with 

£<o. 

• When I > Q > the trajectories present a turning point but f is still positive, thus they 
describe motions of particles on R and anti-particles on L . The center of the hyperbolas 
are now located in the R or L sectors, but are "inside" the hyperbolic trajectory A 2 = 

(j^j) ■ In the R sector, the trajectories [3] with turning point £_ describe particles, which 

enter in R by crossing 7i p with £ > 0, pass through maximum £ at the turning point, 
and quit R by crossing 7i\. The trajectories [4] with turning point £+ connect and 
X^ F , and pass through a minimum £ at the turning. 



When Q < 0, the center of the hyperbolas are located "outside" the limit hyperbola A 2 = 
■ Trajectories [4'] corresponding to the turning points £ + describe also particles 
(because f > 0) connecting Z pp and Tr F - They differ from those with | > Q > by the 
fact that the "partner" trajectories [5] (the other branches of the hyperbolas) intersect the 
R sector. Finally, trajectories associated to the turning points £_ describe anti-particles 
entering R by crossing and leaving by H,\. 

On the quadrants P and F, (e = —1), the sign of the momentum^ = — m£e 2a ^ remains constant 
along the trajectories and allows us to distinguish particle trajectories from anti-particle ones. 
Particles move with £ > on F and £ < on P. At null infinity = ±£oo = qE/2ma > 
in accordance with the fact that particles [6] enter into P from the right past null infinity 2# p 
(with dr > 0, d\ > and £oo < 0) while anti-particles [7] enter from the left past null infinity 
( w hh dr < 0, dX < 0) and £oo > 0). Similarly, on F, particles [8] go asymptotically to 
the right component of the null infinity (2p R ) (with dr and d\ positive) while anti-particles [9] 
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go near Ip L (with dr and d\ negative and d£ > 0). On the horizons, f and Q have opposite 
signs. So, the variable r varies monotonically when Q < but when Q > 0, f vanishes at 
£ = 7^1og[^|r]. When f2 > 0, on the P sector, the corresponding (particle) trajectories [6] 
connect Tp R to Hp and the (anti-particle) trajectories Xp L to Hp. On the F sector, they 
describe anti-particles connecting Tp~ L to Hp or particles between Hp and 1p R . The world 
lines of typical trajectories are depicted on fig. (2). Note that, due to the conformal character 
of this picture, the neighborhood of infinity is contracted into a finite region and the asymptotic 
tangency of trajectories with their asymptotic null rays is no more explicitly apparent, excepted 
in the fact that both reach the same points at infinity. 

3 Field quantization and Schwinger effect 

The quantization of a charged scalar field in an external background electric field is straight- 
forward. It consists into three main steps: 

1) Introduce a covariant derivative = d^ — iqA^ and determine a complete set of solutions 
of the wave equation 

V fl Vy = m 2 (j) , (3.1) 

2) Separate the modes into two classes {0a} associated to particles and {^} associated 
to anti-particle, normalize them using the scalar product built from the usual bilinear 
current 

<— > 

J/i((f>A,(f>A>) = -i(f>* A V V <t>A> (3.2) 

integrated on appropriated Cauchy's surface E: 

(<j> A , <t>A>) = J»(<f>A, <f>A')da? , (3.3) 
and define the quantum field 

$ = J2a(A)<j>A + J2 b+ ( A) r A , (3-4) 

A A 

whose amplitudes satisfy usual commutation relations : 

[a(A), a + (A')] = b AA> [b(A), b+{A')} = 5 AA , 

[a(A),b(A)]=0 etc... 



(3.5) 



3) Build the Fock space from the vacuum state |0) defined by 

a(A)\0) = = 6(i)|0) . (3.6) 
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The interesting physics, the well known Schwinger phenomenon Q], results from the scattering 
of the waves on the external field. In the framework of second quantized field theory, a mode 
describing an incoming particle evolves into a mode associated to the superposition of outgoing 
particles and anti-particles ||. Let us briefly recall how this mechanism works in Minkowskian 
coordinates. In the gauge (|2.2|), equation ( |3.1| ) reads: 



= m> , (3.7) 



but as in the gauge A = Ezdt, the Dalembertian operator and d t commute, the solutions 
of equation (K77p can be expressed as superposition of the modes e~ lat e l ~ tz (p a (z) where the 
functions <p a (z) obey the equation: 

d 2 z + {a + qEz) 2 ] <p a (z) = m 2 <p a (z). (3.8) 

whose general solution can be expressed as a linear combination of parabolic cylinder functions 
Dy[C]; (see for instance refs [|5], Q for a detailed discussion of the solutions of this equation and 
their physical interpretation). The modes describing incoming (anti-)particles are expressed in 
terms of parabolic cylinder functions (whose integral representation will be given in section 6). 

€ m = jjD^ i [C] e~ iCT V^ , ft™ = imL j [-C*] e^V^ 2 (3.9) 

M l 2qE 2 M l 2qE 2 



with ( = e 4 y/2qE(z + ^). Equivalently, modes corresponding to outgoing (anti)-particles 
are given by (j) p a out [t, z] = (<j) p J n [-t, z})* and (j)* out [t, z] = (0^ n [-t, z})* . The constant M has 
been fixed in order that these modes are normalized as follows: 

(€ m ,€n = (<%»*, 4%*) = 5(a - a') , 
(C m ,C' n ) = (€° ut AT t ) = S(a-a') , 

<or*> = (€ mt , €^*) = o • (3.10) 

These two basis are related by the Bogoljubov transformation: 

la out* 



Wl™* - 5 <% m (3.11) 



with: 



V27T nm 2 -n Ttm 2 -tt 

7 = — 5— e 4 9 is e l 2 6 = e 2 * E e 

r[| + <g] 



coefficients that verify the charge conservation relation I7I 2 — |<5| 2 = 1. 
The in and owt vacua are related by: 

|0, Mm&, in) = N~ l/2 e~^ a " b ° |0, Mink, out) (3.12) 
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where we have fixed an arbitrary phase by setting N = (n<7 |t| 2 )- Equation fl3.12| ) show 



explicitly that the m-vacuum is filled with o-ut-particles. This is the content of the Schwinger 
effect. Computing the probability of persistence of the vacuum, we obtain: 



< 0, out\0, in > I 2 



exp 



£rn(l + H : 



e- LTr (3.13) 



i.e., in the limit of large space-time volume LxT, the rate: 

r = — ln(l + e~^£) (3.14) 

2 7T 

of pair creation by unit of space-time volume. The link between the number of modes J2& an d the 
factor L T qE/2 n is obtained as follows || §]. Classically, the trajectories of charged particles 
in a constant electric field are hyperbolas. Each is characterized by the location of its center 
and its radius = m/qE. Quantum mechanically, it appears that it is in the neighborhood of the 
classical center, in a region of space-time extension of the order of (m/qE) x (m/qE) that a wave 
packet describing an incoming particle gives birth to an anti-particle. Accordingly, the number 
of modes relevant to the pair creation mechanism in a given space-time domain, of dimensions 
TxL, are those who centers belong to this domain. The density of modes of frequency a in 
the time-t interval T is (T/2 7r) da. The centers of the relevant modes are located at the 
points z = —a/qE belonging to the space interval L=zi — z 2 - Therefore the values of a, which 
contribute to the sum, are those include in the interval [—qEzi, —qEz<2\ so that in the limit 
of large space-time volume (in order that the boundary effects that are not taken into account 
in this counting of modes become negligible) we obtain J2a = (T/2ir) S„ezi da = T L (qE/2 n). 



4 Rindlerian field quantization 



In Rindler coordinates, the wave equation (|3.1|) becomes: 



em 



(4.1) 



As d T commutes with the Dalembertian operator, the general solution of the wave equation can 
be expressed as a superposition of modes 0^ = (e~ WT / \f2~Tx)T lJJ (^). The functions J-J£) obey 
the second order differential equation: 



em 2 e 2 < 



MO = o 



(4.2) 



S 



The asymptotic behaviors of the solutions of this equation?] are easily obtained by the standard 
W.K.B. technique. The modes 0^0 can be expressed as: 



qE 



1/2 



+ C exp -i (e^ + (u, - ) ( *|Y ' e"< 2_ 

~ V 4a qE*jj\qEJ ^ 



(4.3) 



when £ « +oo i.e. near X + and X , and 

1 



2|w| 



{/^ exp i [u f ] + -D 1 ^ exp -i [w f ] } 



2,71 



(4.4) 



when £ 
by it = 



-oo i.e. near the horizons. More precisely, near X^ and X L , whose points are labeled 
£ coordinates, the relevant part of a mode is given by the weight of its exp(— iuju) 



component, i.e. its coefficient. Similarly, near X^ and X^ , it is its exp(— iojv) part, i.e. its 
coefficient that governs it. Near the horizons, there are the components D°i exp(— iuv) and 
exp(— iuou) that become relevant according to whether we are close to 7i\ or Til on one 
hand or near TCr or Ti^ on the other. Of course all these different weights (C+, Cf_, Df_, D+) 
are not independent. They are related by charge conservation. The scalar product ( |3.2| , |3.3| ) 
is defined on the Rindler patches R and L by choosing as Cauchy's surfaces sections r = C te . 
On such sections the scalar product of the modes we are considering reads as: 



< 



/+OO ri 
(u + u'+^-e 2 ^^^ 
-oo d 

= Qfc L 8{u>-<J) . 



(4.5) 



The prefactor Qr^l is given by the relative sign of the orientation of the vector field d T with 
respect to the global time orientation: 9r = 1,9l = — 1. Using the asymptotic expressions ( |4.3| , 
of the modes we obtain 



Q 



R,L 



&R,L 

8r,l 



|C^| 2 + sgn(^) \D"_\ 2 
\C^\ 2 + sgn{u) \D^\ 2 



(4.6) 



in accordance with the wronskian theorem. By taking as surfaces of integration r = ±oo, we 
easily see that each term of these sums represents a fraction of the total charge of the mode 
that comes from past or goes to future components of the boundaries of the patches R and L. 
More precisely we have: 



Q 



R 



Qu 
1 

Qu 



Q 
Q 



Qi 



n- 



QZ L ' F + Q n t 



Qu 



Q 



(4.7) 



4 The solution of this equation will be expressed in terms of Whittaker's functions in the next section. 

5 In principle, we would have had to indicate, for each mode, the patches on which it is considered, and index 
the various coefficients C and D that it defines by labels R, L, P or F. As we found the notations cumbersome 
enough as they are, we shall omit these precisions, hoping that the context will always be clear enough to avoid 
ambiguities. 
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with the various charges associated to the mode (p^ given by the coefficients appearing in its 

2+ + 

asymptotic expansion: Qu> R ' F = + |C+| 2 , Q h r = + sgn(u;) \D^\ 2 ... 

On the quadrant P and F, the scalar product is given by the integration of the current ( |3.2| ) 
over £ = C te surfaces. It expresses itself in terms of the Wronskian of the solutions J 7 * and 
of eq. (pD as: 

<— > 

< 4>lu,4>lo' >p,f = 8p,f i d $ T w b{uj - uj) 

= Q p j F 8{u-u') (4.8) 

where, due to the time orientation of the vector field we have Op = I, dp = — 1. 
Here the evaluation of the charge is simpler. The Wronskian theorem tells us that: 



Q2* = e P , F sgn(u) 



ID"} 2 - \D1\ 2 







P.F 



ci? - \c u 



(4.9) 



As in eq. (|4.7|) , the total charge Q^' F can be split in terms of the amounts of charge crossing 
the different components of the boundaries of the Rindler patches. We obtain : 



= Qui + Qw = Quf' L + QlT' r 



(4.10) 



with QJ' R = +\C^\ 2 , Q 1 ^ = +|C^| 2 ... 

The particle assignment of these modes is obtained by considering the asymptotic behaviors of 
wave packets of almost fixed value of ui built on these modes, for instance slightly spread around 
a value uj with a Gaussian weight. Indeed, asymptotically, they are expressed as superpositions 
of W.K.B. solutions exp iS(r, £, uj) where S(r, £, uj) is the Maupertuis action for the Hamiltonian 
( p.3[) . By evaluating these superpositions in the saddle point approximation we see that these 
packets are localized on classical trajectories dS/du\- = C te . This allows to interpret them 
in terms of particles, in the light of the analysis of the previous section. Moreover, having 
this particle interpretation of the modes, we will obtain, in the subsequent subsections, a 
precise description of the Schwinger mechanism of pair creation in the Rindler vacua; fig. (3) 
summarizes the results. 



4.1 Quantization on the right quadrant (R) 

In Appendix [A] we discuss the solutions of eq. ( J4.1j ) in terms of Whittaker's functions and, 
using integral representations of them, we give the expression of the various coefficients needed 
to determine the charges carried by the modes. These considerations lead us to define on the 
right quadrant in- and out-modes as follows: 

• Particles in-modes are given by the functions: 

p-iu!T 

V Z7T V 2a 2qE>' l 2a 



2a £ 



(4.11) 
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whose supports (in terms of wave packet) are on Ti R , H R and X R and 



in,R 



in,R) 



27T K 2a 2qE>' l 2a 



(4.12) 



whose supports are on X R) X R and 7^. Obviously Vf n and modes are orthogonal; their 
supports do not overlap on the Cauchy surface H R Ul fi . They are also orthogonal among 
themselves: 



qE 



7TU 

e 2a 



in '^' * a 7 sinh(7rcj/a) 
- a/) 

2 ( W / a 



cosh 



7T 



a 



m 



2qE 



5(lu - J) 



Q(Vr n , R )5(u; - J) . 



—\ e^ + ^cosh 
a / 



7T 



00 

a 



2qE / 



5{uj - uj') 
(4.13) 



The modes Uf nR enter into R via H R while the modes Vf nR emerge from X R . They are 
normalized such that Q(Uf n R ) = sgn(cj) = ±1 and Q(Vf nR ) = +1 represent their total charge. 
From our previous discussion we obtain : 



Qn+( U ?n,R) = sg n M 
In the same way we obtain: 



e _2 fe 


cosh 


" m 2 ' 
7T 2qE 


cosh 




m2 V 
2qE J 



7vm 

e 


|sinh[^]| 


cosh 


H 


a 2qE J 





= sgn(cj) q x 



Q2 



(4.14) 



Qr+(KU) = ?i , Q H+R {V? n , R ) = sgn{u) q 2 

Q{V?n, R ) = QlAVr n , R ) = Qnt^nfi) + Qxi{V?n,R) = +1 



(4.15) 



The interpretation of these equations is obvious. The modes Uf n R describe incoming particles 
when uj > and incoming anti-particles when to < 0. The modes Vf n R always describe incoming 
particles. A schematic drawing of wave packets built with these modes are depicted on fig. (3). 
During their voyage in R, each incoming mode splits into two outgoing branches, one ending 
on X R and the other crossing H R when they leave the quadrant. Charge conservation implies 
that: 

Q = Qh+ + Qi+ 

It II 

Moreover, in a constant electric field only positively (resp. negatively) charged particles can 
arrive from X R (resp. l£) and go up to X R (resp. X R ). This is reflected by the positivity of 
Qx-(Vin) an d Qx + ^in,R}- The charged scalar quantum field operator can be represented on R 

as 



/■+00 

J — oo 



in,R 



(4.16) 
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the various field amplitude operators being defined as usual by: 



<#(w)io)5 



u > 



wio>, 



R 



u < 



«vHlo) 



R 



and the R m-vacuum stateQ as the tensorial product 

io)2" = ff fio)S®io)5 



(4.17) 



Similarly we may consider modes associated to outgoing particles. They are given by: 



V 



out,R 



"■2a 2qE>' 2a 



2a 



(4.18) 



with \N(Vo UtR )\ = \N(U-i n R)\- They correspond to modes entering into R from H R and I R 
and leaving it via 7i R . They carry the asymptotic charges 

QOZ**) = QkM**) = sgnH , Qn-M**) = sgnH q, , Q {V^ R ) = q 2 . 

xi Jx xi 

(4.19) 

These modes are also orthogonal. They describe outgoing particles (u > 0) or anti-particles 
(u < 0) and of course obey the charge conservation rule: 



(4.20) 



with Qj- {V^ ut R ) always positive. As for the m-modes, to obtain a complete set on R, we have 
to add to the set of functions ( |4.18| ) the modes 



out,R 



V Z7T H 2a 2qE>' L 2a 



(4.21) 



which are also orthogonal among themselves and with respect to the set {V^ ut R }. These modes 
come from 7i R and I R and end on I R and X R . They charge content is 



These charge are such that: 



out) 



(4.22) 



(4.23) 



With these modes we may express the quantum field operator ( 4.16| ) as 



+00 



+e(-u)b+™ t (u;)V^ R } . 



(4.24) 



6 This vacuum state is not the end of the story. Indeed the true vacuum state must be defined on a Cauchy 
surface of the whole space, i.e. it must include a factor representing the L part of the vacuum. 
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Of course out- and m-modes are related by a Bogoljubov transformation. Using formulas 
recalled in Appendix we easily obtain : 



u: 



out.R 



in,R 



V 



out,R 



(4.25) 

9{u) {a&(w)Z^ + a£vMKU} + 9{-u) {tvvHK,* + e^ u {u)U^ R } (4.26) 



whose coefficients, fully displayed in Appendix B, are such that: 



R 
UV 



:<o)\ l 



-vu 



\Q 



R 
VU 



KHI 2 = I^MI 2 = l7^vMI 2 

and (remember that qi and q2 are function of u): 



q-2 



(4.27) 



o&(w<0)| a -|/^( W <0)| 1 



a* v (uj>0)\ 2 + \a« u (u J > 



R 



qi - q2 

qi + <?2 



(4.28) 



The last relation shows that the Bogoljubov transformation between the in- and out- vacuum 
state factors is non trivial only for u < 0, in which case we obtain 



io>5 io)5 



<*wvM e Q «v(") " ^ >\0) \Q) (Vw<0), 



(4.29) 



reflecting the instability of the in-vacuum and the pair production predicted by the Schwinger 
mechanism. The fact that only the u < vacuum factors appears to be unstable is in accord 
with formulas ( 4.14j) ,( [4.15| ) which indicate that it is only for these values of u that an in-mode 
of a given charge generate owt-branches with opposite charges. Note also that this instability 
can be interpreted as a tunneling between classical trajectory living in R and its "partner" 
when this latter visits the R quadrant. 

The probability of vacuum persistence is given by: 



R R 

I < 0,in\0,out >" K I 2 



J| \a U v{u)\ 2 

-£i*(i + fel 2 ) 



w<0 

exp 



L uj<0 



exp 



uj<0 



_ m 

1 + e""W 



q E e ^ 



(4.30) 



The significance of a symbol as XL<o is analyzed in ref. ||, § and recalled at the end of 
section 3. The density of modes of frequency u, in an interval of width T in the r variable, 
is T/27T did. The sum over u < corresponds to consider essentially the particle trajectories 
whose mirror trajectories also visit the R quadrant. The spacelike coordinate £ c of the center 
of contributing trajectories are obtained immediately from equation ( |2.9|) which gives du = 



^^dQ = —aqEp c dp c and thus T d\u\ = qEdV. As a consequence we obtain: 



| R < 0, in\0, out > R | 2 = exp 



In 1 



m4\ qE 



qE 



2ty 



V + T x finite terms 



(4.3i; 
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where in the limit of large space-time volume V (i.e. large interval [— U\, —u)2\ of length Auj ^> 
a ), we have 



V 



r 

qE' 



-Au 



T 

2a 



(exp[2a£ c (-w 2 )] - exp[2a^ c (-a; 1 )]) = T R 



(4.32) 



= Pc(^2) — Pci(^i) and introduced 
^From eq. ( 4.31 ) we see that the 



where we have denoted the difference of Rindler radii by R 
a mean proper time T = T(p c [u2) + p c i(^i)]/(2 a) = T R. 
rate of Rindler pair creation from Rindler right quadrant vacuum consists in two pieces. The 
first one, as expected, is proportional to the measure of the volume multiplied by the Schwinger 
rate. To this dominant (in the large volume limit) term, a surface correction, which can be 
expressed in term of dilogarithm (Spence) function Li 2 , has to be subtracted: 



£1+6" 

u<0 



T 

2^ 

aT 

Air 2 



L0 2 



1 + e "Tfie" 22 ^ j duo 



As emphasized by Robert Brout [10|, it is instructive to compare the formula for the per- 
sistence of Rindler vacuum, eq. Q4.30|) , with a similar situation which arises in the theory of 
Hawking black hole radiation. In particular, our result stands in strong analogy to the analysis 
of s-wave emission by a Schwarzschild black hole, and even more so when this latter is reduced 
to a pseudo (1 + 1) dimensional problem, when one neglects the effects of the s-wave "cen- 
trifugal" barrier (See ref. 0]). It will be noted that the rate eq. ( |4.30|) defines is a difference 
between a volume term, proportional to R, and a surface term wherein R is replaced by a 
characteristic length a/ (qE R). Each of these terms is proportional to T, hence giving rise to a 
rate per unit volume and rate per unit surface respectively. For the Schwinger case, the object 
of our present study, R is macroscopic and the surface term is negligible. But for the case of 
s-wave black hole evaporation, the steady state production comes from a region of 0(M) of 
the horizon where 2 M is the Schwarzschild radius. Recall the formula for the rate of particle 
emission 



(T r ) 



71 



, , -k B T^ wk . . (4.33) 

J- Hawk. ±^ 

There is a prefactor M in the above, the volume of the "skin"; this replaces R in our formula. 

Equation ( f4.33p contains the vacuum expectation value of the energy flux in Unruh vacuum, 
the outgoing mode piece of the Hartle-Hawking vacuum. In ref. [O there is displayed the 



same calculation in Boulware vacuum, the analog of what we have called Rindler vacuum. 
It is shown that in this case vacuum polarization effects cancel against the emission effects 
and <T;> 

Boulware 



0. In our eq. ( [4.31 



this subtraction comes up in analogous fashion, 
due to the particular structure of the Rindler modes, and hence of Rindler vacuum, in the 
vicinity of the horizons. The difference between the two cases is that in the black hole case the 
cancellation of (T or ) Boulware must be total. This is a consequence of the static character of the 
Schwarzschild metric which underlies the construction of Boulware vacuum. In the Schwinger 
case the problem is actually intrinsically time dependent. Indeed the electric field has to be 
(adiabatically) switched on and off in order to prepare the m-vacuum initial state (or to observe 
the crat-one). 
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4.2 Quantization on the left quadrant (L) 



Formally the field equations and their solutions are the same on the left and right quadrants. 
Their only difference results that on L the vector field d T points near the past. This implies 
that the signs of the charges and the in and out labels have to be interchange with respect to 
their values on the R quadrant. The analytical expression of the in- and out-modes on the 
L-quadrant can be deduced immediately from the drawing (see fig. (3)) of the supports of 
the modes built on the R-quadrant. Adding a extra subscript (R, L, etc . . . ) to make the 
distinction between the modes defined on the various quadrants, we obtain for the m-mode the 
expression: 



-IU)T L 



Uin,L{iL,T L ) 

OC Kut, R {ZL,T L 

and for their charge content: 

Q-H+ i U in,L 
Ql+ (Min,L 

Similarly 



H2«t 2qE>' 2a 



2a 2 6 



(4.34) 



—Q T - (ptout,R) 



- sgn(w) q 2 
-Qi 



Here again, we may define in- and out- vacua, Fock spaces, etc 
persistence is given by an expression similar to ( 4.31|) . 



(4.35) 

(4.36) 
(4.37) 

, and the amplitude of vacuum 



4.3 Quantization on the past quadrant (P) 

On the past quadrant, it is the vector field — d^ p that defines the future direction. The various 
modes can again be expressed again in terms of Whittaker's functions: with their charges given 
by eqs ( |4.8|) . More precisely we may choose as particle incoming modes 



V? n AtPiT P )=Af(V 



in,PJ 



-lUJTp 



'2ir 



W 



V 2a 2qE > ' 2a 



qE 



—i — -e 

2a 2 



2a 5p 



They came out from I R , carrying a charge 

Qx R {vt n , P ) = W{vt n , P )\ 2 

and leave P via Ti^ and 7i£, taking away the charges 



/ qE\ _z_(<±-n£ 
I- — Je 2(a i E 



+1 



a 2qE J 



cosh 


IT — 


sinh 


TTUJ 

a 





sga{u) q 3 



(4.38) 



(4.39) 



(4.40) 
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and 



Q W7 (vr„,p) = -e-^ 



cosh 


L (»- 

7T - 

_ \ a 




m 2 y 

2qE)_ 


sinh 


■KLO 

a 





= - sgn(u) q 4 



The complex conjugates of the anti-particle incoming modes are given by 



Z€,P(£P,TP) = A/Wn,P 



-lUJTp 



oc 



. 2tt 
V^p^p, -tp) 



*(:£t 2qE^' l 2a 



2a2 C 



Their asymptotic charge content, when they appears on X L , is: 

Qt-(ZC,p) = -1 • 

This charge splits into 



(4.41) 

(4.42) 
(4.43) 

(4.44) 



(4.45) 



with sgn(u;) (q 3 — g 4 ) = 1. 

Outgoing W-modes on P have to vanish on Hi i.e. their coefficient £>_ must be zero. They are 
given by Whittaker's M-functions: 



lUlTp 

e" a?p M 2 

271 ( 2a 2qE>' t 2a 



'2a 2 



with 



iAr(^,p)r 2 = i^ + (-Mi)i 2 = M 



e 2a 



(4.46) 
(4.47) 



whose charge content is 



Q(Kut, P ) = Qn R (Kut,p) = +mnn 
Qr-RV) = -94 <0 , 
Qr-RV) =93 >0 . 



(4.48) 



Accordingly, modes with positive (resp. negative) values of uj have to be associated to particles 
(resp. antiparticles) and, as expected, the sum of the two last charges reproduces the total 
charge of the mode. Finally, we have also to consider the modes 



-tWTp 



V^, p (£p,tp) = m^ ut , P )^=^e-^M ^ 



2qE>' "2a 



qE 



2a£_ P 



OC 



(4.49) 
(4.50) 



that emerges from X and leaves P via H R , with the charges 

Q(Kut, P ) = Q wz (V D V) = -sgnH 

Qx- L (Kut,p) = "93 , 
Qx-(Kut,p) = 94 • 



(4.51) 
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Here also the quantum field operator can be expressed in two equivalent forms, using the in- 
and out- modes: 



+00 

—00 
+00 



did 



duo 



0{u) «(w)ZCp + b+^V^p 



out,P 



(4.52) 



(4.53) 



Using the relations given in Appendix [S], we obtain easily the Bogoljubov transformation 
relating in- and out-modes: 



whose coefficients (see Appendix B) satisfy charge conservation relations: 



and: 



K>>o)f-ltt>o)f 



huvi^ < 0)| 2 - \e£ u {uj < 



-1 



(4.55) 
(4.56) 



From these equations we may repeat the steps leading to the probability (|4.31| ) of no pair 
creation and obtain: 



,P , P ,9 

< 0,in\0,out > 



exp 



exp 



Eln(l + I/^| 2 ) + 5>(1 + 



-v p I 2 
luvl 



Loj>0 



LU<0 



^ , „™ 2 cosh7r(^ 



m ' 
2qE- 



uj>0 



sinh 



7TCJ 

a 



1 + e 



2qB y 



cosh 




sinh 





(4.57) 



4.4 Quantization on the future quadrant (F) 

The modes on the F quadrant are obtained from those defined on P in the same way we pass 
from R modes to L modes : 



out,F 



(£f,Tf) oc Ut p {£ F) t f ) 



(4.58 



the proportionality factors reflecting the arbitrary phases appearing in the normalization fac- 
tors. Note also that, as expected, these modes satisfy charge conjugation relations 



ZCt.F^F, Tp) OC V^ p(( F , -Tp) 



(4.59) 
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Figure 3 summarizes all the discussion of this section. We have schematically represented, for 
the four Rindler quadrants, the behavior of typical wavepackets built with the different modes 
discussed here above and the charge content that each of them carries when it crosses the 
horizon or reaches infinity. 

5 Vacuum decay rates revisited 

In this section we shall reobtain the rates of vacuum decay by using functional methods. Our 
motivation is that we hope in this way to obtain an alternative picture of the physics underlying 
the pair creation. In addition, as such methods often necessitate the use of approximations, 
it is useful to use them as a test in examples for which the exact answer is known since there 
sometimes arise physical situations in which they are the only tool available. First we shall 
briefly recall the general formalism and illustrate it in the context of our scalar field interact- 
ing with a constant electric field, the calculation being performed with respect to an inertial 
t, z coordinate system, i.e. with respect to the vacuum states defined by inertial observers. 
Here, the interaction being quadratic, all the integrations are gaussian and the calculations 
can be performed analytically. Then we shall consider the same problem with respect to ac- 
celerated observers, living in the R quadrant. We shall first perform our calculations using 
the expression of the Feynman propagator as sum of modes and recover the results obtained 
in the previous sections. Then we shall consider the Schwinger propagator which expres itself 
as a kernel integrated with respect to a time variable (Schwinger proper time 7 ). Thanks to 
the invariance with respect to translations both in time t and r (boost) we shall be able to re- 
duce the expression of the Schwinger propagator to the evaluation of path integrals of ordinary 
one-dimensional quantum mechanical problems. We shall then evaluate these path integrals at 
one-loop approximation (which is exact for quadratic potentials) and compare the results to 
those previously obtained. Finally we shall also compute the in-out vacuum amplitude, using 
an expression for the Feynman propagator obtained by unwinding the inertial propagator, and 
show that it does coincide with the amplitude calculated in the secure framework of the mode 
analysis. 

The previous vacuum persistence amplitudes obtained from "mode" calculations, may also 
be derived from the Green functions of the quantized field. Indeed, as is well known (see for 
instance the book ||IT[), the vacuum persistence amplitude can be expressed as a functional 
integral : 

(0, out\0, in) J=J , =0 = Z [0, 0] = e lW (5.1) 
where Z[J, J*} is defined as : 

Z [J, J*} = J V(f)V(j)* exp iiS [0, 0*] + i J J*(f)d d x + i J J(j)*d d x^ (5.2) 

in terms of the action S [(f), 0*] of the charged scalar field, minimally coupled to the electric 
field: 

S [0,0*] = ~ J ^{(D M 0)(^0)* + (^0)*(^0) + (m 2 -^e)00* + (m 2 + *e)0*0)} 
7 Also called fifth time in the framework of quantum field theory in 3+1 dimension. 
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dxdx'(p*(x)H xx '4>(x') . (5.3) 

(Here, and in the following, we denote collectively by x, x', the coordinate variables of the field; 

the context indicating if they are Rindlerian or minkowskian.) 

A standard computation allows to give a sense to the expression of W : 

/*oo P 

W = -ztrln(-£) = / dm 2 d d xG F (x,x) (5.4) 

where Q = —TiT 1 is related to the Feynman propagator 

G F {x,y) = -i(0, out\T {cj)(x)(i) ] \y))\0,in) / (0, out\0,in) (5.5) 



by: 



g(x,y) = ±{G F (x,y) + G F (y,x)} • (5.6) 



Thus, the imaginary part (the one which encodes the vacuum instability) of W reads: 

yoo p poo r 

ZmW = Ke dm 2 / d d xG{x,x) = Ke dm 2 / d d x G F (x, x) . (5.7) 

Our purpose now is to evaluate this expression for an inertial observer and an accelerated one, 
using different schemes of calculation. 

5.1 Mode representation of the Feynman propagator 

If the asymptotic vacua are those of an inertial observer, using (t, z) coordinates, it is straight- 



forward to obtain the Feynman propagator [12] as superposition of modes (eqs p.9| and followings 
in the text) : 

iG F (x,x') = - [ daip p a out (x)<p a a out (x') 



7 

+6{t - £') / da^ out {x) V p a out *{x') + 6{t' - t) I da V a a out *{x) V a a out {x') (5.8) 



Therefore we get : 



G(x,x') = -^-^^(^(^^(xO-^*^)^*^)) 

f da (<f%»*{z)tfr**tf) - rt«**(x)tir*tf)) . (5.9) 



2 

The first term in this last expression is traceless, so we obtain : 

5_ 

7 



Tie tr Q =lm- J d 2 x J da <p p a out (x)(p a a out (x) . (5.10) 



19 



• Using the integral representation a la Schwinger for the product of wave functions derived 
in the Appendix C, we can rewrite this as: 

Tie tr Q = -l m e —^ f £ x f da f°° *i -im^E^ft^Es (5 n) 

2tt 2ct J J Jo ( S inh2g£s)2 

and the imaginary part of the W amplitude reads: 

1 y/qE f 2 f f 00 d.S e ~im 2 s ,„ TrJ ^_ s _^ i 



ImW = Im-^ / d 2 x / da / , "* — e W+%)^«Es ({U2) 

2tt 2iir J J Jo (s-tfa) (smh2qEs)? 

After the t and z integrations we get, in accord with eq. (|3.13|) : 

T r . r°°- i7T d6 _ j2 £ T 



ImW = -— da — — e~ l ^ y = — ln(l + e""^) . (5.13) 
87r J J-oo-in 9 smh 9 An 

• In Rindler coordinates (r, p) = (r, a^ 1 e a ^) on R, the in-out Feynman propagator, ex- 
pressed as a mode superposition, reads as : 

iG*(x,x') = 9(t-t')M™ du{U™\x)U™ t \x')+V™ t (x)V™ t *(x')} + 

r ° du} ^u: u \x)vT{x') 

-co a.ux> 

+ 0{r'-r)f cLu^-U^(x)V^*(x') (5.14) 

J-oo Otuv 

and the part which encodes the pairs production is given by : 

ftetr£(x,x) = 1m f d 2 x f° duj—U™\x)V™*(x) . (5.15) 

J J-oo &UV 

Using the integral representation (given in Appendix C) of the product of modes appearing in 
this last equation we get, once the space-time volume element is reexpressed as d 2 x = apdpdr: 



f f° 1 ™ 2 ™ f°° 
Tie tr Q(x,x) = Im dr du—e 2 i E e 2 °- / pdpe 

J J-oo An Jo 



50 ^ e 'C*-*V*^V ia (-^) . (5.16) 
oo cosh u a 2 cosh u 

This representation allows us to perform the integration on m? and we get for the imaginary 
part of the functional W: 



ImW = ^JLxm [drf du [°°^ -JL-e^'eW [°° pdpe^e^^e^ 
An J J-oo J-oo-iZ 6 smh 9 Jo 



a 2 sinh 9 



where we have introduced the new variable 9 = u — i~ and rewrite an Bessel-7 function as a 
J function. The next step consists now to carry out the p 2 integration which, thanks to the 

leads us to : 



formula (6.611.1) of [|T3 



ImW 



Im ■ 



Air 



dr 



o 



du 



d9 



sgn(Ke6) ^ 
9 sinh 9 



l e i ^ e {sgn(ne9) sinh 9 + cosh 9} 1 
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Finally, by splitting the 9 integral according to the sign of Tie 9 we obtain: 

2 2 

ZmW = — dr — I duo — — dr — I du- 



o7T J J — oo 9 J-oo sinh# 8 71 J J — oo 9 J-oo sinh# 

= — / rfaj ln(l + e"~) - — / rfcu ln(l + e"~e 27r ~) (5.18) 

47T J-oo 47T J-oo 

in accord with our result ( |4.30|) established by using the Bogoljubov transformation coefficient 
between in- and owt-modes. 



5.2 Schwinger representation of the Feynman propagator 

For the sake of completeness, let us briefly recall the so-called Schwinger representation of 
the Feynman propagator. Formally, the Feynman propagator Gf is the inverse of the kinetic 
operator K = 5 jS^Ml {—T> U T>^ + (m 2 — it)} and can thus be computed as Gp = —K^ 1 = 

—i J °° dse~ tKs . To express it, Schwinger Q introduces the kernel defined as K(x,y;s) = 
(x'\e~ lKs \x) which obeys a Schrodinger equation : 

(V^~m 2 + t-^-)K(x,y;s) =0 (5.19) 

with the boundary condition K(x, y; s — > + ) = 5 d (x—y). Using this proper time representation 
of the Feynman propagator 

rod 

Gf(x,x')= / dsK(x,x';s) , (5.20) 

JO 

the vacuum persistence amplitude can be written as: 

/*oo p /*oo 

lmW = lm dm 2 d d x dsK(x,x;s) . (5.21) 

Jm 2 J JO 

The kernel K(x, x';s), defined as the matrix elements of the evolution operator K, can be 
expressed path integral: 

K(x, x'; s) = J vX»(s')e iS( - x ' x '-' x ^ s '^ (5.22) 

where the domain of integration covers all the paths that connect the point x to the point x' in 
a (rescaled) time^ s, and S(x, x'\ x^(s'); s) is the classical action along these trajectories, given 
by: 

S(x,x';X^(s');s) = ds\^- + qX»A^X[s\) -m 2 | . (5.23) 

Moreover, let us emphasized that if it is often supposed that the Feynman propagator obtained 
from modes coincides with the Schwinger representation ( |5.20|) displayed here above, it is only 
in a few cases that their equivalence has been proved (See for example |L4|, |i5|). 



8 This time variable is proportional to the proper time measured along the natural trajectory that connects 
x to x', but, in general, does not coincide with it. 



21 



Reduction of functional integrals 

The Schwinger propagator introduced above via a functional integral is a formal object that 
has to be defined more accurately. To this end we remind the Feynman prescription |[{|. We 
divide the time interval s into iV subintervals of length e, such that s = Ne, and introduce N 
times the closure relation in the expression of K(x,x'; s) taken as the matrix element: 

K(x,x';s) = (t',z'\e- iAs \t,z) 

= J (j[d Zi dt^ (x'\e- iHe \x N ) . . . (x x \e- i&e \x) 

r ( N \ N 

= / J] dz i dt i II - tj) 6 ( z j+i - z j) - ie (xj+i\H\xj}} , (5.24) 

J \i=l I j=0 

before taking the limit iV i— > oo at the end of the calculation. 

Let us first consider the problem for an inertial observer. Here the hamiltonian operator 
reads: 

H= {(d t -iqA t ) 2 -d 2 z +m 2 } = [-(p t - q A t ) 2 +p 2 z + m 2 }. (5.25) 

We compute its matrix elements in position representation, with \x) = \z) (g> \t). Using the 
orthonormalized eigenvectors of the momentum operators p z \k) = k\k) and pt\oj) = —u\u>) 
which connect to the position eigenvector by (k\z) = e ^- z and (t\u) = e J^ - Thus, 



(x j+ i\H\xj) = -5(z j+ i - zj) I dw(t j+1 \u)(u\(pt - qA t (zj)) 2 \tj) 

+ 5(t j+1 - tj ) Jdk(z J+1 \k)(k\p 2 z \z,) + 5(z J+1 - Zj)8(t j+1 - t. W 
= -5(z j+1 - Zj ) J ^(w + qA t {z 3 )) 2 e-^-^ 



+ 5(t j+1 - tj) J ^-k 2 e ik ^-^ + 5{z 3+l - Zj)5(t j+1 - tj)m 2 

e ik(zj + i—zj) e —iu(tj + i—tj) 

dkduj H d (k,u,Zj) (5.26) 

Z7T Z7T 

where we have introduced the classical Routh function H c i(k,u, Zj) = —(u + qA t (zj)) 2 + k 2 +m 2 . 
Inserting this in eq. (|5.24 ), the kernel K(x, x'\ s) reads as: 

N N ,. ikj(zj +1 -Zj) p -iu)j(t j+1 -tj) 



K{x, x';s) = {H dzjdtj) H / dkj duj ^eH^k^^) 

r N rlh r rli > r N 

= j(J[—l-) j ri° e -two(*'-*) J Y[_dzjt ik ^e~ itHcl( - k ^^ hZi) 



j=0 j=l 

N "du. 



8=1 

duo 



Y[ dZi j ^ e --o(*'-t )(4me) -^ e ^/ Ve( . +(? ^) 2e -W 
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where the last equality defines the Fourier transform of the kernel K(x, X J s\ ctS the path integral 

K u (z,z?;s) = f Vz(s)e l fo c ^ ds \ (5.28) 



with respect to the (classical) Lagrangian C u (z, z) = i 2 /4 + (uj + qEz) 2 — m 2 . 

A similar calculation can be done in Rindler coordinates (r, p = a" 1 exp(a£)), where we use 
the t independence of the hamiltonian to perform the reduction of the path integral to a one 
dimension quantum mechanical problem. However, due to the explicit coordinate dependence 
of the metric, the calculation do not reduce just to a strict rewriting of the previous one. The 
Schwinger kernel is now given by : 

K(x,x';s) = (p',r'\e-^ s \p,r) 

J ( HdpidrA J] {KPo+i ~ Pj)K r j+i ~ r j) ~ ie(x j+1 \(ap)*H(ap)*\xj)} 



\fd 2 pp' 



The appearance of extra pa factors results from the coordinate dependence of the volume 
element. Here the hamiltonian reads as : 

«=P?-(p,-^ 2 ) 2 ^-^ + ™ 2 (5-29) 

where the momentum operators are p T = jd T and p p = \{d p + ^-). Compared to the classical 
hamiltonian, an extra quantum term l/4p 2 appears, just as for the plane rotator |fl5| |. This 
is dictated both by the requirement of hermiticity with respect to the measure p dp and the 
classical commutation relations (see refs |17]])- The orthonormalized eigenstates of these op- 
erators and the connectors between them and the position eigenvectors (\x) = \p) <8> \t) , with 
(p'\p) = {apY 1 5{p ! - p)) are: 

p p \k) = k\k) , p T \u) = -uj\u) , 

(k\p) = (ap)-i^= , ( r | w ) = . (5.30) 

V27T V27T 

So we obtain the matrix elements: 
(x j+ i\(ap)^H(ap)^\xj) = (a 2 pjp j+1 )^ (x j+l \H\xj) 

$(Pj+i - Pj)( T j+i\(Pr - qEA T ^)) 2 \Tj) 



a 2 p 2 



r e ik (Pj + l~Pj) I 

+ <H r i+i - T j) J dk — k 2 + (m 2 - -^)S(r j+ i - T j )S(p j+1 - pj) 

e ik(p j+l -pj) g-iwfa+i-T,-) 



4pf 



dkduj H c i(k,u},pj) (5.31) 

Z7T Z7T 

where we have set H c i(k,u, p) = j— -^^{oj + aqE 2 p ) 2 + k 2 — + m 2 . Inserting, the matrix 
elements ( |5.31| ) into the expression of the kernel ( p.29| ) we obtain : 

i r N r N e ik i(pj+i-Pj) e -*^j( T i+i- T j) 
K(x,x';s) = (a 2 pp')~2 J J| dpidn J JJ dkjduj — — e" 



(2 l\ 

[a pp) 



i=l " j=0 



J 1 J 2tT „ 27T 



i=l j=0 
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As in the calculation done in minkowskian coordinates, the independence of the hamiltonian 
with respect to the time variable makes its conjugate momentum conserved. Each r 3 - integration 
generates a delta function of uij — cJj_i, making all uij integration trivial but the first one on 
t^o- Moreover the integrals on the kj are gaussian, and we easily obtain: 

(a 2 pp')2 J 2n i= i j=0 



(a 2 pp')^ J 27r 

where, as for the inertial propagator, we have introduced a Fourier transform with respect to 
the Rindler r time as the path integral : 

K™ n (p,p';s) = J Vpexpz jT £ w (p, p)ds' ; (5.34) 

the lagrangian being now given by : 

1 E 1 

C^(p,p) = p 2 /4 + —(uj + q-ap 2 ) 2 + — -m 2 . (5.35) 
a z p z 2 4p z 

To pursue our analysis we have to evaluate the two (unidimensional) functional integrals 



( |5.28| ) and (|5.34| ). The first one can be evaluated explicitly, but we did not succeed in expressing 
the second in a close unambiguous form. To go ahead, we shall evaluate them at one-loop 
approximation, which actually is exact for gaussian integrals and thus gives the correct answer 
for the functional integral referring to minkowskian coordinates ( |5.28| ). In this approximation 
we just have (in principle) to compute the classical action S c i, as the well known Pauli-Van 



Vleck formula |Tj| [TP| gives : 

K(q, t; q', t') = J Vq(t) exp iS(q(t)) = det g^) h exp iS d . (5.36) 

In principle this formula necessitates the evaluation of the hessian matrix d 2 q ,S c i. At first sight 
this implies the knowledge of the general two point expression of the action. But actually it 
suffices to know the general classical solution of the equation of motion to obtain it. Indeed it 
can be shown (see for instance ||25|| ) that the hessian matrix is the inverse of the matrix built on 
Jacobi fields vanishing at point q. These are simply obtained by varying the general solution 
of the equation of motion with respect to the initial velocity componentsP|||T9|. 27 



Let us apply this method to compute the kernels (|5.28| and |5.34|) relevant to the evaluation of 
trace of the Green function and rate of vacuum decays. Here we have only to take into account 
closed paths in the functional integrals, starting and ending from the same value of the z or 
the p coordinate. 

In z-coordinate the calculation is straightforward. The classical motion equation is: 

£_ V *(, + ^) = (5.37) 



9 The reader will easily recognize in this method the n-dimensional generalization of the method advocated 
by S. Coleman to compute functional determinant pqj. 
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whose particular solution such that z(s' = =F§) = z is given by : 

Z[S) ~ qE + [Z+ qE } coshqEs ' 
The value of the classical action computed along this closed trajectory reads: 

t s 

3d = I ds' £ u (z,z)\z=z(s>) 



u 2 cosh2qEs - 1 2 
= qE{Z+ ^ sin^qEs ~ m3 " ^ 

Here the computation of Jacobi fields is immediate, because for quadratic potential it only 
consists in solving a harmonic oscillator problem : 

' -Aq 2 E 2 )h(s') = . (5.40) 



\d 2 s 

whose solution obeying the Cauchy conditions h(s' = — |) = and i^?h(s')\ s > = -2. = 2 (because 
the unusual normalization 1/4 of the kinetic part of the lagrangian) is: 

h(s') = -^-smh2qE(s' + J) . (5.41) 
qE 2 



Thus we get: 



K w (z,z;s) = ( . . ^ s " lh2 ' £s . (5.42) 

27U smh 2g£/S 



Inserting this expression in eq. ( |5.21| ), the vacuum amplitude persistence, after t and m 2 
integration (T — J dt), amount to : 

l.qE iT r f°° ds _■ 2 „ f dz iaE sMiqEa (z\ " ) 2 

-f- ^2 / ///,) / P %m S / cosh 9 £ S <,E J 



ImW = lm-{ — )2— / rfcj / _ e - jms / j-e^ssariB;^^ 

% 2m 2tt J Jo s J (smh.2qEs)2 

= — I dwlnfl + e - *^) . (5.43) 

47T J 

This one loop computation gives the exact result, as expected, because here the lagrangian is 
quadratic and the one-loop approximation (|5.42|) is actually exact. 

Let us turn now to the same problem formulated in Rindler coordinates, for the accelerated 
observer. We have to evaluate the functional integral ( 5.34 ) on closed paths starting from and 
ending on p(s' = — ~) = p(s' — |) = p. The equation of motion, derived from the lagrangian 
is: 

1 d 2 p 2 .uj 2 1. 1 , n2 n . AA , 

2d£ + + 3 f " 2 q E P = ° <5 ' 44) 

which, by use of the energy theorem, reduces to 



2 



dp\ 2 Q 



2 



q 2 E 2 \ds' J ^ vr p 2 
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2(3+(p 2 + — ) (5.45) 




where (3 is an integration constant and 

"=^ + i >0 • (5 - 46) 

the correspondence between this quantum number and the classical constant of motion being 
Q fa (m/qE) 2 Q. When (3 > 0, obviously there is no turning points in p; the motion goes 
from p = to p = oo, and the trajectory never crosses twice the same point. Actually, to get 
a trajectory with a bounce on the potential barrier, (3 must be less than —Cl. Hereafter we 
restrict ourselves to this class of trajectories. Integrating equation ( |5.45| ), we obtain as solution 
of equation Q5.44|) : 



pV) + (3 = e(pl + fl ««h{2g^- «>) + ¥>} (5.47) 

cosh if 

in terms of arbitrary constants po, (3 and sq and where e = ±1 and we have set, for further 



notational convenience, cosh tp = (p§ + (3)/\j (3 2 — Vt 2 . The solution which bounces in time s, 
such that p(s' = — | ) = p(s' — | ) = p is: 



f W) + P ={p ^ + pfJ^Ml (5.48) 

cosh<p 

where the value of (3 is now fixed to f3 = jp 2 — cosh(qEs)^J p 4 + Q 2 sinh 2 (gi?s)j / sinh 2 (g£'s), 

which implies that <p = —(qEs). Obviously (3 is negative and p 2 + (3 will be negative or positive 
according to p 2 is less or greater than fl i.e. according to the position of p with respect to 
where the minimum of the potential p 2 + Cl 2 / p 2 occurs. The computation of the value of 
the classical action along this trajectory is straightforward. We obtain, in terms of the variable 
if), defined through the relation sinh^ = f2sinh qEs/ p 2 : 

, VuJ 2 + a 2 ( cosh qEs - cosh iff \ 
S u (p,s) = r-7— , (5.49) 

on which one it is easy to check that: 

esL^l = _tm (5 , 50) 

It remains now to compute the Jacobi field h(s') obeying the Cauchy conditions h(s' = — |) = 
and d s ih(s')\ s i = ^s = 2. Varying the general solution (|5.47|) of the equation of motion with re- 
spect to the "energy" parameter (3 we get a Jacobi field vanishing for s' = — |. Once normalized 
such that d' s h(s') f s = 2, its value for s' = s/2 becomes: 

s 2 

h(-) = — ^sinh qEs cosh-?/; (5.51) 

and gives the value of the Van Vleck determinant. Collecting all these results, we obtain for 
the path integral Q5.34D with p = p', the approximated expression : 



2 i ■ qEuj i ■ V dj^+a^ / cosh qEs — cosh lb , ( \ 

— tm Z S-\-l^ S+l— . r ; -+W) 

p a a \ sinn ip T / 

Kf n {p, p; s) « = , (5.52) 

' -^sinh qEs cosh if> 
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and for the vacuum persistence amplitude 



_ _ T r , r°° ds .-a.,^, qEVL r°° I( (coshV>) 2 , v^w f ^B- co ^ ,\ 

ImW = Im du — e im 3+1 a s \ l / dip- -^e ° ^ sinh ' /, ' 

2ni J Jo s V 16z7r Jo (sinh'0) 5 

(5.53) 

Let us evaluate the if) integral by the saddle point method. Indeed, given the approximations 
used to evaluate K^ m (p,p; s), it would be inconsistent to endeavor to go beyond this approxi- 
mation. The phase reaches its maximum (with respect to if>) at if) — qEs, corresponding to a 
saddle point of width (Vto 2 + a? /a) cotYiqEs. This leads to the gaussian approximation : 



Jo r (sinh ip)^ \ u 2 + a 2 sinh qEs 



, (cosh if>) 2 ■ y^ 2 +a 2 / coshfl-cosh^ ,\ I 2ina e* v qEs 
dip— -Ve - ^ sinh * + ^ « W— — (5.54) 



and the vacuum persistence amplitude reads: 

ImW ^ Tm— { du T e " im2 V^V"^ 9 ^ (5.55) 

47ri J Jo s sinh g£s 

« Jm-^ / du I" r^— f^(cu)e- im2s+i2 ^ s + ^(-u;)e- im2s ) . (5.56) 
47u 7 Jo s sinh g£s ^ / 

Comparing eq. ( |5.56|) with eq. (|5.18|) we see that their difference reduces to 

Rc /^Lw fc -2? (5 - 57) 

once we take into account the i e prescription. This divergent quantity is physically meaningless. 
It is independent of the physical parameters of the problem and must be subtracted in order to 
recover the stability of the Rindler vacuum in the absence of the electric field. Of course all the 
finite corrections we have mentioned before, and that we recover here above, are mathematically 
meaningless; the approximations we used for the Schwinger kernel making them all irrelevant. 
Our motivation for having discussed them is just to show that no spurious term are introduced 
in the dominant contribution to Im W. Actually this leading order can be obtained directly 
by ignoring all subtleties but just using a quadratic approximation of the potential : 



p z a 2 4p z 

- m 2 ) + q 2 E 2 (p - p ) 2 (5.5* 



around the minimum p§ = fl m -^oj. The computation is similar to the one discussed for the 

inertial observer once we have substituted m 2 — — to m 2 . As for the evaluation of (|5.42j ) 
we obtain now : 



whose trace gives expression ( 5.55 ), analog to (|5.43 ). 
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5.2.1 Is the zero winding propagator the Rindler propagator? 



In the appendix of ref. |J, the Fourier transform of the inertial Feynman propagator was 
expressed as a sum of winding terms, similar to those given in ref. ||28| . Each term of this 
sum corresponds the Fourier transform of a path integration restricted to the subclass of paths 
joining their ends after having performed a fixed number of winding around the common vertex 
of the four Rindler quadrants (around the origin O, see fig. (1)). It was conjectured that the 
kernel K™ =0 obtained from the zero winding sector, provides the kernel leading to the Rindler 
propagator : K™ =0 = K Rm . Hereafter we reinforce this conjecture by showing that the rate 
calculated from this zero winding kernel : 



K 



w=0 



1 e 



-i(m 2 -^)Es 



An sinh qEs 



, q e pp' 

2 sinh qEs 



(5.60) 



leads to our previous results. To this end we substitute in eq. (|5.21|) the expression of 
K Rm (x,x, s) obtained by using the expression ( |5.60| ) in eq. ( |5.33| ). Reexpressing the Bessel-7 
function as a sum of Bessel-f^ functions and integrating over the p variable, we obtain : 



ne^ f °° p dp Ho du / °° kt°(p, p; <0 f 

= TZe£ JZo du / °° e-^s+^es)^ f oo p dp 



1 1 E „2 



qEpe ~T 



p coth qEs 



2 7r sinh(^_E s) 



gg p 



2 sinh(<j_E s) 



+ (2 0(cu)sinh(7r^) - e^) (- 



' 2 sinh(<j_E s) 



)} 



(5.61) 

Formula (6.611.3) of ref. JTB| provides us the result of the p integration. So we obtain for the 
vacuum persistence amplitude, computed using fl5.60| ): 



T 

XmW = TZe— 



du 



An J-oo sinh(7r^) Jo s smh(qE s 



ds 



9{uS) sinh(7r— ) e 
a 



%2^-qEs 

n. * 



10 , 



— 9{— lu) sinh(7r— 
(5.62) 

i.e. expression ( |5.56| ) which is equivalent to eq. ( [5.181 ). 

Let us note that the Fourier transform of the Schwinger kernel has to satisfy the equation 



' A p A 



i 2+ A 



K^p,p';s) = 



(5.63) 



\imK w (p,p'; s) 

s— >0 



(5.64) 



and the initial condition on the R quadrant 

5{p - p') . 

The expression ( |5.60| ) satisfies this equation and the initial condition in the following sense: 

\p,p';s) =5(p-p , ) + [2^Msinh(7r-)-e^]5(p + p / ) . (5.65) 



lim^ =0/ 



This extra term is not a surprise. It is localized on the boundary of the R patch and reflects 
the correlations between the left and right quadrants that appear when the vacuum decay and 
that pair creation occurs. Actually the existence on the horizon of a singular contribution to 
the Rindler (uncharged) scalar field has been put into evidence some time ago by Parentani 
29| who has shown that such terms are essential to the validity of the theorem stating that the 



usual Poincare invariant vacuum state is the state of minimal energy. 
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6 Construction of Unruh modes 



In this section we shall establish connections between Minkowskian and Rindlerian modes. 
Instead of working with integral transforms, we find more convenient to diagonalize the Bo- 
go ljubov transformation. This is achieved by introducing Unruh modes. They correspond to 
new basis of the Fock spaces, which share the same quantum numbers as the Rindler modes 
but that do not mix positive and negative frequencies of the Minkowskian modes. The pur- 
pose of this section is to build these modes which allows to easily describes algebraically the 
Minkowskian particle content of the Rindler vacua. Completeness of Minkowskian and Rindler 
sets of modes insures each element of one set can be expressed as superposition of members 
of the other one. To obtain the relation between them we shall, instead of computing directly 
the overlapping integrals, made use of integral representations of the parabolic cylinder func- 
tions (see ref. ]13 , eq. (9.241.2) for the mathematics and [g] for their physical meaning) . 



Accordingly, we may express the in- and out- modes (eqs |3l^ and followings in the text) as : 



^ 1 ) r[§ - V5F(W /4 



o 



dxe ixV ^ E( ^- + ^~ t ^)x~' 1 ^-^ (6.1) 



1 e 4 e 4< ? B 

roc 



2 



€ m (U,V) = —. 1— r- — 



r dxe-^^ + ^-^)x~ 1 ^ (6.2) 
Jo 

where we have used the ingoing and outgoing light-like coordinates U = t — z and V = t + z. 
Using the reflexion operation U t— > —V and V \— > —U, we obtain : 



>C ut {u,v) = (<f a m (-v,-u)y 
% mt (u,v) = {<f) a a m {-v,-u)y 



On the other hand, from the integral representations ( A.5|) of the Whittaker's function, we get 
the two representations of the Rindler modes (see eqs. |A.12j ): 



w:(u,v) 



( qE i3 E UVf qEV\ (Jt e-fe^ 
' 1 e * ( ) 2« — 



Uvra 2 / 



' 2qE ' 2 

m 2 1 



qE \ 1/2 el ^ UV( qEU\^ 



Ana" 



feZ. -e — 



L y l a 2qE ^ 2' 



dte iU r uvt t i -- % ^-^(i + ty^ 

(W+(U,V) 



m 1_ 

2 



(6.3) 



(6.4) 



(6.5) 
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where we have reexpressed the Rindler coordinates £ and r in terms of light-like Minkowski 

coordinates, through the relations £ = ^ ln(— ea 2 UV) and r = ^ ln(— e]0. 

Let us now consider a wave packet, denoted hereafter as fi^ m , of inertial in-particles 0£ m : 

S>',. ■ / ^(-p^.-W-'JC-^) . (6.6) 



30 (J Cl> (J 



Using the integral representation ( |6.1| ), we can rewrite this superposition as the convolution : 



e im e f e -^v(v-u) 



JxyfZqE 



V-U :X Z 



2 x 2 i E 2 F 



10 

a 



qE(-V + X\ 



qE' 



dx 



(6.7) 



where we introduced a two variable function, F, built from the coefficients of the wavepacket 
(|6.6D through the relation: 



a 



/+0O i.) -2 

A(a, -yP° e - ia - da 
-oo (X 



(61 



If we choose the arbitrary function F such as : 

F 



a 



\2a 2 J 



(6.9) 



(the choice of the coefficient in front of this expression anticipating a subsequent normalization), 
eq.( |6.7| ) becomes an integral representation of Rindlerian V m-modes on P and L , i.e. W~ 
functions (see eq. (| A . 1 2|) ) : 



2„2 \ 4 



.in 



57T a 
~qE 



ese s i E 
M 



(6.10) 



While general arguments of completeness insure the existence of the function A(a/y/qE,u/a), 
we need its explicit expression to obtain the continuation of on the rest of space-time (in 
the regions V > 0). Taking the inverse Fourier transform (with respect to the variable (3) of 
eq. (|6.9|), we obtain : 



g 2qE POO 



/r qE J a 



/qE_y 
\2a 2 J v^Wo 



e" 2 e V9 E p '« + '2?E 2 dp 



(6.11) 



r[i( 



m 



zgil/ a 2 p) 5 yqE 



3j7T 0" 



V2aV ^ 

Using twice the integral representations ( |6.11| ) of the parabolic cylinder functions, the in- 
rals over and x become straig 
functions, i.e. modes lA in ^ and M ut,F 



tegrals over -y=j an d £ become straightforward. We obtain an expression involving only W, 



qE 
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(6.12) 



More explicitly, by introducing the characteristic functions Xl,r,p,f (which are equal to zero 
or to one according to that the coordinates considered correspond to a point belonging to the 
sector labeling the function or not), we may write : 



2 /cosh[<-^, 



—i xf e 2a 



m a;' 

2gE a ■ 



sinh(^ 

ll 



_m_ _i_ J.1 

2qE a ' *~ 2 J 



rr 1 _ a m 

1 h L 2qE\ 



in,R 



+ 



rr 1 — o m 

1 h l 2qEl 



out,F 



(6.13) 



The mode (|6.13|) so constructed defines an Unruh mode. Indeed it is a (superposition of) 
positive frequency modes, algebraically related to Rindler modes. Similarly, we may define 
other orthonormalized m-Unruh modes through wave superpositions built as in eq. ( |6.6| ) : 



— 



•pm 



a 


10 




a 


a 


uJ 


VqE' 


a 



a 



iqE' 



e 1 ^ sgnM 



(6.14) 



a 



XpU?* P -ie 



sinh(^; 



rU(22Z- — —\ -L- 1 

XnUin, R + « 2 ° — ^ r^ai ^ V ™,l 



r[±-< 



2g£J 



pr-/ mi _ on _,_ ll 
1 1-2 2qEi 



3 # , 



XLUr* L + iXFe^Ur* F 



(6.15) 
(6.16) 



where the phase ip is given by : 

■0 = arg 



r[^l 



r[^(S + i) + J1 



(6.17) 



Let us emphasize some properties of these Unruh modes. By considering on I~ the supports 
of the Rindler modes that define them, it is obvious that these Unruh modes are orthogonal to 
each other and normalized. This also results from the relation : 



a 



d(^=)A(a,u)A*(ea,uj') = 6(e)8{u'-u) 



iqE 



±1 



(6.18) 



that can be easily proved using the integral form ( |6.11| ) of the coefficients A(a, u). On the 
other hand, once we realize that these modes have to correspond to Minkowskian particle 
or antiparticle m-modes, starting from I R or they may be built from purely algebraic 
considerations. Indeed it suffices to fix their behavior on I~ and continue them across the 
horizon by requiring continuity (always in term of wave packets). For instance, suppose we 
want to start with modes that born on Ip R , we have to take on P , Vi n ,p modes (see fig. (3)). 
Then in order to maintain zero Cauchy data on Xp P and X~[ we have to paste them respectively 
to Uin^R and Vi n ^ modes. Then it remains to fix the linear combination of modes defined on F 
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that fit our construction on Ti^ and Tl^. So, using the relation ( |A.7| ) we may fix the relative 



weights (and phases) between these modes on each quadrant. And an a priori unexpected 
property appears: the interferences are such that the particle Unruh mode Q pm vanishes on 
Tp L and the antiparticle Unruh mode Q ain on Tp R . This finds its mathematical origin in an 
analyticity property: the fact that the w Unruh modes can be expressed as linear combination 
of only functions Aif and the modes Q as pure combination of Wf functions. 

Just as we have introduced an in Unruh basis ( 6.6|6.16j ) we may consider a out Unruh 
basis built out from superposition of out-particle and antiparticle modes with as coefficient the 
complex conjugate of those used in the construction of the in Unruh basis. Actually these modes 
are obtained by performing a space-time inversion [(tf,£f) ( t p,£,p), ( t r,£,r) (tl,£l)] on 
the previous modes and taking their complex conjugate. They are given by : 



n 



p,out 



pout 



a out 



a out 



-oo 
oo 



a 



qE a 

° U \*Jj>out 



a 



qE 
a 



iqE a 
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d(-TTS) 



qE' 
a 



iqE 1 a 



d( -=) 



qE' 



a 



U \*iaout - a 



'qE a 



-)*cj) a a out d(- 



iqE' 



(6.19) 
(6.20) 
(6.21) 
(6.22) 



These modes are related to Rindler owt-modes by a Bogoljubov transformation given in Ap- 
pendix B and schematically represented on fig. (4). 

To compute the coefficient of the Bogoljubov transformation between the in and out Unruh 
modes, we use the Bogoljubov transformation ( |3.11| ) to transfer the relation between modes 
into relation between connectors linking Rindlerian and Minkowskian modes. For instance we 
obtain for the modes QP° ut : 



Of. 



out 



x \/qE ^qE 



(6.23) 



Fortunately, the 7 and 5 coefficients are a independent and factorize out of the integral. The 
contribution of the antiparticle modes sums directly (see eq. ( |6.16j )) and gives : 



- 7T 7rm 



(6.24) 



The summation of the particle modes needs a little bit more work. Thanks to the linear relations 
connecting parabolic cylinder functions, we get: 



. 7T _ 7T771 7TCi 



prl _|_ i<a — 

A(±a, uj)* = — ° 2gE \ e~^e~^e^A(To-, u) + e^e~^A(±a, uj) 



2,71 



(6.25) 



and the particle modes contribution to eq. (|6.23| ) can be written as a sum of two in Unruh 
modes. So we get: 



out 



_ 7rm ttuj 



iy/EV + e-^e^^qle^'w p J n + ie~^w a w 



(6.26) 
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where q\ is defined by eq. ( |4.14| ) and ip' is the phase: 



In the same way, we obtain: 



ip' = arg 



I LH a 2qE 



T[\+i 




(6.27) 
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(6.28) 



We found instructive to recover, using these basis, the decay rate (|3.13|) between in and out 
Minkowskian vacua. Formally we obtain : 



< 0, owtlO, in > 



exp 



T f°° / _irm 2 

2 — / du In 1 + e " E 

27T J-oo V 



(6.29) 



Let us emphasize that 2{T/2n) /f^ da; = 2(X^>o + X^<o) gives the expected total space-time 
volume factor LT/2tc, in accordance with our interpretation of XL>o an d XL<o as the space- 
time volume of the various quadrant. 

More interesting is the computation of Rindlerian particle population in Minkowski vacua. 
Expressing the field as superposition of Unruh modes and Rindler we obtain, thanks to eqs 
( |6.13| - |6~T6| ), we obtain the Bogoljubov transformation between Unruh and Rindler creation and 
annihilation operators : 
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(6.30) 



the various coefficient a, (3, 7 etc . . .being defined implicitly by eqs ( |6.6| , |6.14| - |6.16l) . These 



relations allows us to express [|30| the in Minkowski vacuum as (see Appendix B): 
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(6.31) 



and can then be interpreted as a superposition of Rindler's pairs consisting in Rindlerian par- 
ticles and antiparticles and their partners, anti particles and particles localized in the opposite 
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sector. Using the explicit form of the Bogoljubov coefficients given in Appendix B, we obtain 
the mean density number of Rindlerian out particles that an accelerated observer (in the sector 
R) detects when the quantum state is the in Minkowski vacuum : 

= (0,Mink, inlaid* (uj)a™*(uj)\0, Mink in) 

-7rm 2 

= e~~ (6.32) 



and: 



n K ut R = ( Q , Mink , in\a\^\u)a^(uj)\0,Minkin 



(1 + e " E ) 

e a, — i 

for the various particles, while the density number of antiparticles is : 

n V out oR = (0, Mink, in\b\;™\uj)b ^(uj)\Q , Mink,, in) 



(6.33) 
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(6.34) 
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In the same way an accelerated observer in the L sector will detect the following populations: 

n V out = riifout = e <* E (6.36) 



of antiparticles and: 



_ 

(1 + e ^e~ 2 "-) 

n K%L ~ n vr<'o« ~ e -2^ _ i • v°- 6 <) 

of particles. 

All these populations are given by Bose factor and corrective terms proportional to the Schwinger 
factor. On the other hand, an accelerated observer in the right sector measures a total charge 
given by: 

roc /'OO r-0 

Qr — / duN V out + / duNyout — / duN V out 

JO wR J-oo bjR J-oo wR 

= / due~^ . (6.38) 
Jo 

Charge conservation implies that a left observer sees: 

r°° m 2 

Ql = - I duoe~^ . (6.39) 
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Note (a check of all this algebra) that, always due to charge conservation, the same expressions 
of Qr and Qp are obtained if we compute them using in density number of particles, but the 
individual contributions of U and V type of (anti)particles are different. 
Of course the total charges Qp and Qp are zero. 

Moreover, when — > (i.e. for a weak electric field or massive particles), some populations 
( |6.32| , |6.36| ) vanish, and the other ( |6.33| , |6734| , |6.35| , |6.37|) become Boltzmanian in character m 

l/e~^~ — 1. This can be understood as follows. Let us consider, for instance, the U° u p modes. 
In the limit where qE goes to zero, these modes become localized near i° R . Their charges 
(see eqs (|4.22|) ) on the horizon components vanishes, while on the past and future horizon 
they become equal to unity. Indeed these modes tend to Bessel-I functions, i.e. function 
which grows exponentially when p goes to infinity, and that do not contribute to the Hilbert 
space of the neutral quantum field (see for instance ref. Classically, they correspond to 

hyperbolic trajectories pushed away to infinity, hyperbolic trajectories whose radius (m/qE) 
become infinite and whose asymptotic points on and Zj£ slip to i°. On the contrary, V° u p 
correspond classically to hyperbolic motion whose past and future asymptotic points tend 
respectively near i~ and i + i.e. inertial trajectories. When such trajectories cross the R 
quadrant, they enter in and leave it out across the horizon components Tip and Tip]. This is 
confirmed by the fact that for such modes their charge contents (eqs ( |4.19| )) vanishes at infinity 
but become equal to unity on the horizon components. 



7 Uniformly accelerated charged detector 

In a previous work M we have computed transition amplitudes (see hereafter, formulas [7.26| to 
|7.34| ) between charged particles of masses and charges that we denote here M, Q and M', Q' 
(instead of M, Q and m, q as in ref.0) interacting by exchange of a third kind of particles of 
mass m and charge q (that was called fi and a in ref.0). Their interactions were described by 
a three field interaction Hamiltonian and first order perturbation calculations were performed. 
Our purpose now is first to clarify some aspect of the physics of the accelerated detector and 
to make an explicit contact between it and the model built on the three interacting fields. 
The detector can be seen as a "two-level ion" propagating in a constant electric field. The ion 
levels are supposed to have rest masses M and M' (resp. charges Q and Q'), both much greater 
than the mass m (resp. charge q) of the exchanged quanta. This ion makes transition, without 
recoil, between its two levels, i.e. whatever is its internal configuration, it always moves with 
constant acceleration a — Q E/M = Q' E/M'. Hereafter we show how this model appears as 
a limiting case of the three field model, by comparing amplitude of transition obtained in both 
cases. In the following, all the calculations are done at first order of the interaction coupling 
constant. The physics of the accelerated two level detector is described through the effective 
hamiltonian : 

H int (r) = ag [Ae~ tAMr e iq I A ^ dx " ^ (t,^ = 0) + A^e iAMT e~ iq I A ^ }(r, £ = 0)] , (7.1) 

where we have taken into account the prescribed trajectory of the detector £ = 0, (p = a -1 ). 
The operators A and A^ are annihilation and creation operator acting on the two dimensional 
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Hilbert space generated by the two quantum states up and down of the ion 



A\up) = \down) , A\down) = , A^\down) = \up) , A^\up) = 

(7.2) 

The energy difference between these two levels is given by the mass gap AM > 0. Indeed, for 
the uniformly accelerated detector following the trajectory £ = 0, the local rest frame time ( 
the detector proper time ) coincides with the Rindler r time, fixing the conjugate energy to 
be the rest frame energy difference between the two levels. The integral phase factor / A^ dx^ 1 
insures gauge invariance of the interaction. It has to be computed along the trajectory of the 
detector, from an arbitrary origin to the position of the detector at interaction proper time r. 
For the model we consider here, it reduces to up to an arbitrary constant of integration. 

Before discussing this model, we would like to recall the reader a few results concerning 
the simpler situation where the exchanged agent is massless and chargeless. In this case the 
field can be expressed, on R, as a superposition of Unruh U and V modes (see ref. [|J for 
details) : 

0i? = 4> U r + <I> V r = J o du \a% <p% + aZvZ) + herm. conj. (7.3) 

The detector transitions, when the field state is the vacuum, can only occur by emission of 
Minkowskian neutral quanta. At first order of perturbation, the probability amplitudes of such 
transitions from the lower mass state of the detector to the higher one (denoted by B) or the 
converse (denoted by A) via the emission of one quantum associated to an Unruh U mode are 
given by : 



B(u; AM) = (0 M \(+\a w %e~ l S dTWnt \-)\Q M ) = ~i~9^J(^ + 
A(lu;AM) = (O M \n auj T T e- l I dTWnt \+)\0 M ) = -ig J^ a J(*f - 



(7.4) 



where a w = 1/ — e -22 ^) and T T denotes the time ordering operator with respect to the 
detector proper time. Its is interesting to compare these formula with their analog where the 
emitted quanta are those associated to Minkowski modes of energy k (eq. (2.48) of ref. ||). In 
the latter case the instant of emission is given by f(k) = a~ l lnk/AM with a width 1/yAM. 
When the emitted quanta are in pure Unruh's states, their "energy" are fixed (u = —AM), 
but the process becomes completely delocalized in time. In other words the Rindler "energy" 
balance is actually a "boost" constant of motion conservation relation. When recoils effects 
are taken into account but the classical picture of the detector still valid, the ion jumps from 
one hyperbola to another whose centers are such that their difference in localization remains in 
accord with eq. (|2.9| ) for the different values of u. Summing the square of these amplitudes and 
interpreting as usual the 5(0), we obtain the detector excitation probability per unit of proper 
time : 

PB(AM) = J du\B(u-AM)\ 2 = ? * M ^ t l — (7.5) 
and similarly for the desexcitation probability per unit of proper time : 

PA(AM)=^ (i _; 2 ^ ) (7.6) 
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The ratio of these probabilities is: 



PB(AM) ——am 

— e a 



(7.7) 



PA(AM) 

which reflects the thermal equilibrium of the detector and the radiation, at Unruh temperature 
a/2ir. 

Now let us repeat the same calculation, but in the framework of a charged exchanged 
agent. First we compute the probability amplitude of spontaneous excitation of the detector 
by emission of an Q (resp w) antiparticle of quantum number to. Expressing the quantum field 
as a superposition of Unruh modes, we obtain at first order in the coupling constant : 



B 



; AM) 



-iga(0,Mink,mrt\bp> (uo) / n mt (r)dr\0, Mink, in) 
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£j-ja out * 



g a 5 
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NqE 7 



J dre^ AM -^ 



1 „out*r^ ~\^out , ,\„pout 



Here N =< 0, Mink, out\0, Mink, in > is the normalization factor (|3.12j ), 5/7 is defined by 
eq; fl3.12| ) and 



out 



and e°r u n are the coefficients of the (non frequency mixing) Bogoljubov 

\ zd r 

transformation between Unruh and Minkowski modes ( |6.19||6.22"D . Finally, expressing the Unruh 
modes in terms of Rindler modes (eqs P-12|) , the r integration in eq. ( |7.9|) becomes trivial and 
gives 5 functions. So we obtain the excitation amplitude probability B(fl u ; AM) of the detector 
by emission of an f2 w quantum as the sum of 



— %- 

and the integral : 
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(7.10) 
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(7.11) 



Using the relations ( |6.18 ) and (|6.25|) for integrals of products of A{a, uo) functions, it easy to 
show that the a integral gives a S(uo — Co) function; as a consequence the integral over Co is trivial 
to carry out. So we get for the second term in the expression (7TS) of AM) : 

gy/2n6.,AM qE uo^A uo m 2 . 1 

—i o{ - — )r( — 1 — 1- 1 — 1 — 
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2a 2 



(7.12) 
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which, once the explicit expressions of the coefficients Cqv an< ^ a ?w R ( e 1 s are explicited, 

lead to: 



, 37TOJ 7TT7T 
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(7.13) 



In the same way we obtain the excitation amplitude via emission of a vo^ antiparticle : 
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(7.14) 



which, again using the relations (|6.18 , |6.25|) , lead to : 
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Similar calculations lead the desexcitation amplitude of probability of the detector, by 
emission of a particle of type Q or w : 
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(7.16) 



and: 



A n7 u ; AM = -*^r^£ fea + 7 V " : 
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On the other hand let us notice that the transition amplitudes due to absorption of quanta : 

B'(^ ; AM) = (+|(0, Mink, out\Te~ l f drnint ^a^ (u)\0, Mink, m)\-) 

1 • 
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A'( QuJ ;AM) = (-|(0, Mink, aa^Te^ I '^w^M&jpt (a?) |0, Minfe, in) |+) 



are directly related to the previous ones by : 

B'( flu ;AM)=A( flu ;AM) and A^" ; AM) = ; AM) 

tUu W M W,,, •57,,) 



(7.17) 
(7.18) 

(7.19) 



These last relation can be verified by repeating the previous calculation, or more directly using 
CPT transformation arguments ||. 

The occurrence of the Dirac distribution 5(^- — |4 + -) in all these amplitudes makes trivial 
the estimation of probabilities of transition per unit of proper time. With obvious notation, T 
beeing a large detector proper time interval, we obtain : 
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(7.21) 
(7.22) 
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with cu* = AM — qE /2a fixed by the mass shell condition dictated by the delta functions and 



2qE 



where we have set \C\ 2 = g 2 ^§2 ( ^f)~ l e T ^ E / (cosli7r(^ - cosh 

In the limit of small charge exchange ( ^| i— » 0), approximating, in the limit {qE \— > 0), 
Whittaker functions by Bessel functions (see ref. P2"l , section 13.3) : 
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it is easy to see that the probabilities PA(Q, AM) and PB(Q, AM) are of order one compared to 
the probabilities PA(w, AM) and PB(w,AM) which of order e _7rm l qE . Indeed an amplitude 
such as B{wu\ AM) corresponds to the probability amplitude of excitation of the detector 
by absorption of a quantum described by a mode, i.e. a rindlerian mode 14™* on the R 
quadrant. But, as discussed at the end of the previous section (see eq. (|6.32|) ), such modes 
corresponds to fluctuations produced by the Schwinger mechanism, and in the limit of vanishing 
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charge, their populations become exponentially dampen (here we suppose that qE/m 2 « 1) 
compared to the populations of fluctuations due to the Unruh mechanism. Accordingly, the 
probabilities of transition involving such modes are also exponentially small, compared to those 
refereeing to process involving the modes created by the Unruh mechanism. This illustrate how 
the two kinds of creation process cooperate to the physics of the detector. This situation is 
similar to the physics of the Reissner-Nordstrom black-holes which start to loose their charge 
by the Schwinger mechanism, and after by Hawking evaporation [03 1. In the limit of small 



charges we obtain as ratio of the probabilities of transition the thermodynamical equilibrium 
expression : 

PB(AM) = [PB(Q; AM) + PB(zu; AM)] ^ ^ {AM _^) . . 

PA(AM) [PA(tt; AM) + PA(zu; AM)] ~ 1 ' ' 

whose meaning is discussed at length in On fig. (5) we have plot this ratio of total 
probabilities divided by the Boltzmanian factor exp[— (27r(AM/a — qe/2a 2 ) = exp[— 27ra>*]. 
This figure illustrate the role played by the Schwinger factor which control the instability of 
the quantum state of exchanged quanta. For large mass or, equivalently small charge, we 
recover the usual Boltzman equilibrium formula, otherwise new equilibrium relations have to 
be considered. Once again the same physics occurs when we consider a quantum state in a 
strong gravitational field. For instance if a black hole is enclosed in a box, it may reach an 
equilibrium configuration with surrounding radiation. However the radiation alone will not 
obey a local equation of state linking energy density and pressure, but it is only the whole 
system that will satisfy an equilibrium relation involving its total mass energy, the geometry of 
the box, etc. 

In ref [[/J we also obtained expressions for the transition amplitudes of an accelerated ion, 
coupled to a massive and charged scalar field $ m of mass m and charge q = Q — Q' , when its 
recoil is taken into account. The two level ion was modeled by introducing two scalar fields ^ m 
and of masses M and M' and charges Q and Q' . The transitions were dictated through 
the interaction hamiltonian : 

H int = gjdx (& M y M ,$ m + * M tf M ,& m ) ■ (7.25) 
We get for the desexcitation amplitude : 

A(k\k',k") = (0,out | cC(k") a$(k')TeWf^ nt a$(k) | Q,in) 

= C8(k- k' - k")A{k\k',k") (7.26) 

which at first order of a perturbation expansion reduces to 

- c i f k f ' 2 h" 2 \ ^2 
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2 q a 



which X denoting the complicated expression: 
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and where we have set: 
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while the constant C is related to the normalization factors and Bogoljubov coefficients through : 

-ig 1 



C — 

Similarly, the amplitude of spontaneous excitation for the two level ion, defined as: 
B(k'\k, —k") = (0,out | cC{k") bZK-k) Te- i S dtHUU bS(-k') \ Q,in) 

was obtained into the form: 

B(k'\k, -k") = C8(k-k' -k")B(k'\k,-k") 



(7.3i; 
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with at first order: 
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(7.34) 



These amplitudes were obtained by using modes solving the filed equation by separation 
of variables in the gauge A = —Etdz, modes such that their spatial dependence are of the 
form e lkz . To make contact between these amplitudes and those obtain in the framework of the 
detector model we first have to change our Fock space basis in order to use, in both situations, 
the same quantum numbers to label the states. These k modes, gauge transformed so as to 
solve eq. (|3.7f ), are related to the modes Minkowski modes introduced in eq. (|3.9| ) by the 
(trivial) Bogoljubov transformation: 
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(7.35) 
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In the following, we only discuss the first order expansion of the A(k\k', k") amplitude : 
A(k\k',k") = -i(0,out | aZ\k") a°$(k') fdtH int a!^{k) \ 0,in) 



-i I dada'da"a{a" , -k")a{a' , -k')a(a, k)(0,out \ <C\a") a^(a') / dtH int a!^{a) \ 0,in) 



and show how from this amplitude we may recover the transition amplitude of the two level 
detector. 

We shall evaluate the limit of this "3-field" amplitudes of transition for M, M' — > oo with 
AM/M = (M — M')/M -C 1 and Am/ a <C 1, which corresponds to the situation where the 
three field model mimics the heavy two level detector, without recoil. The detector states are 
given by sharply localized states built out of the vacua \0m > and \0m> >• They are obtained 
as : 



+ >= 



J f + (k)dka\ 1 (k)\0 M > \0 M > > , h >= / f-{k')dk'<J M ,{k')\U M > \0' M > (7.37) 



and the operator A\ acting on this two dimensional state space, is: 

A ] = f f + (k)a j M {k) dk f f_{k')a} M ,{k')dk' (7.38) 



These states are supposed to be normed (<+|+>= / / + /* dk = 1, < — \— >= J /_ /* dk' = 
1, < +|— >= 0) and, as in the large mass limit, the Schwinger process vanishes we have not 
to distinguish between in or out vacua \0m > and \0m' >■ At first order of perturbation, the 
amplitude of transition between these states becomes 

A(k") = J f+{k)dk J f.{k')dk'A{k\k',k") (7.39) 
~ -ig < 0,m, out\a™Xk") dt dz <fi M (t, z)$* m (t, z)<p* M ,(t, z)\0,m,in > (7.40) 



where TLintir) is given by eq. (|7.1[) and where <pM{t,z) and <f>M'{t,z) are classical solution 
of the wave equation, given by the superposition of modes weighted by the function f±(k); 
for example: <f>M(t,z) = J f+(k)(jrf.(t, z) dk. Note that at this first order of the perturbation 
expansion, a necessary condition for A(k") to be non zero is that the wave packets 4>M{t,z) 
and 4>M'(t,z) overlap, but at higher order virtual detector states allow "tunneling" transitions 
between non overlapping configurations of the detector. The heavy detector corresponds to the 
limit where these wave packets have their supports blend into a single classical trajectory, here 
a hyperbola whose center is at the origin of the (t, z) coordinates. Such configurations can be 
obtained by using, in the large mass limit, the maximally localized packets introduce in |J or 
more directly just by considering suitable superpositions of W.K.B. approximate solutions of 
the wave equation : 

«M) * (f^*^-«V^ . (7.41) 
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In the limit considered, we thus obtain: 



£ iAM r e -i 



(7.42) 



because AM a = (Q — Q')E = qE. So we recover, from the field description, the energy 
difference between the two detector levels as predicted by the "equivalence" principle. The 
amplitude reduces to : 



A{k") 



i / d</'a(</\-tf')(-\(0,m,(mt\a% t (o'') / dr^ nt (r)|0, m, in)\+) ,(7.43) 



and using the Bogoljubov transformation (|B.12|) between Unruh and Rindler modes, we get : 



A{k") 



e i E 
^qE 



{a^(a 7 u)A(n^ AM) + a^{a,u)A{w w - AM)} 



da 



du 



(7.44) 



which, can be explicitly integrate (thanks to the integral representations of Whittaker and 
parabolic cylinder functions) to give at the end: 
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e 2a 



(7.45) 



where // = A M - Q'E/2a. 

This result can also directly be obtained from the expression ( |7.27| ) of the amplitude 
A(k\k', k") by taking its limit for infinite masses and charges. Indeed, in this limit, the ampli- 
tude factorizes into a term involving only the k and k' momenta and a term that tends to eq. 
( [7.451 ), the effective coupling constant g being then implicitly defined as the product of g times 
a double integral involving the functions f+(k) and f-(k') and the k, k' term that comes out 
from the full amplitude. 

First let us discuss the factor I (eq. |7.28| ). The limits of the parabolic functions D are simply 
obtained by direct substitution: Q \— > : Jf= and £ \— > — i + ■ For the hyper geometric 
functions occurring in I, we use its series expansion and obtain : 



1 1 qE , 

hm 2 F 1 {ie M +-,ie m +-, l+ie M -ie M '\ —p^TB, 

M,m^oo 2 2 Q hi 



-is®, .qE._i-o. 



iqE s 



2a 2qE< 2a 2CL 



(7.46) 
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while for the second, we first have to use the inversion relation to pass from the argument Q' j q 
to q/Q' and the Stirling formula to finally get: 



11 Q'E, 

hm 2 Fi{ie M + r, ze m + -, 1 + ie M - ie m ] ^r, 
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(7.47) 



T{\-it m + i^Y qE ' v ~2a 2 ' 

Moreover this can further simplified using eqs (|A.22 ) and the limit of the A amplitude reads: 
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(7.48) 



This is in perfect agreement with eq. ( 7.45|) if the coupling constants g and g are connected by a 
(q, m) -independent relation, dependent of the precise form of the weight functions characterizing 
the detector. 

So we may relate the terms of the amplitude of transition (|7.27 ) to Unruh modes, and we 
see that the physical reason that favor one kind of modes with respect to the other found its 
root in that their own existence is a manifestation of the Schwinger mechanism and that their 
probabilities of interaction is weighted by a Schwinger factor. 



Conclusion 



In this paper, we have studied, in Rindler coordinates, the quantization of a charged field 
interacting with a constant electric field. The main characteristic of this problem is that it 
involves two acceleration parameters: the acceleration of the Rindler observer (the detector of 
section 7) and the natural acceleration of the charged quanta (qE/m). So we obtain a toy (but 
exactly solvable) model for the quantization of a charged field in a Reissner-Nordtrom black 
hole geometry, in the same way as the Unruh detector mimics the physics around a Schwarschild 
black hole. It is that similarity that constitutes the main motivation of our work. 
Now, let us summarize the main points of our analysis: 
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The quantization of a charged field in Rindler coordinates illustrates the "symmetry 
breaking" between particles and antiparticles in the R and L quadrants. Eqs ( |6.38 ) and 



( |6.39D show that the Minkowskian vacuum state carries a positive charge on the R and (of 



course) the opposite charge on the L quadrant. This was expected, as antiparticles (resp. 
particles) are always obliged to leave the R (resp. L) quadrant, which is not necessarily 
the case for particles (resp. antiparticles) that may stay inside for ever. Nevertheless, 
although obvious, we find interesting to see how this property emerges from the rules 
of quantum mechanics and is encoded in the wave function of the various modes used. 
Let us note that when the quantization is performed on the full Minkowski space, no 
such charge polarization effect comes into evidence because in this framework operator 
expectation values are obtained by averaging over the complete space instead of over just 
one quadrant. 

The discussion of the classical trajectories that leads to interpret the conserved quantum 
number uj as the invariant distance A from the center of the hyperbolic trajectories to 
the common vertex of the four Rindler quadrants shows that (once more in terms of wave 
packets), it is near the center of the trajectories that the pair production mechanism 
occurs. Indeed A depends crucially on the sign of u, the extra term 1/2 being only the 
reflect of the quantum indeterminacy of the position. For instance, it is only for uj < 
(eqs ( |4.30| , |4.31|) that the Bogoljubov transformation is non trivial (i.e. mixes particles 
and antiparticles) on the R quadrant. 

On the Rindler quadrants R or L the persistence of the Rindler vacua is given by the 
usual Schwinger result modified by a surface term ( |4.33|) . The latter becomes negligible 



in the case of constant electric field in the large volume limit but it plays nevertheless an 
important role. In the framework of black hole physics it becomes the vacuum polarization 
term that cancels the Hawking radiation flux in the Boulware vacuum. 

The limitations of the W.K.B. evaluation of the Feynman propagator are exemplified by 
evaluating it on the one hand as a mode superposition and on the other hand using the 
Pauli-Van Vleck approximation of the Schwinger kernel. Here the Jacobi fields were an 
essential simplifying ingredient of the calculation, which shows that such approximation 
is nothing else than a quadratic expansion of the classical potential around its minimum 
(eqsp38j p39|). 



The interpretation of the usual Minkowskian propagator as a sum over winding Rindlerian 
propagators is reinforced. We have shown that the rate of particle production calculated 
from the zero term of the Minkowskian propagator winding number expansion coincides 
with the rate obtained from the Rindler propagator. 

But on the other hand the calculation of the Rindlerian population of the Minkowski 
vacuum leads to mode densities that are not equilibrium distributions in character. This 
is not surprising as in a constant electric field there are no physical reasons to expect an 
equilibrium distribution. 

Moreover some of the distributions ( |6.32| , |06| ) obtained vanish when the electric field 
goes to zero. This reflects the fact that for a charged field we obtain twice the number of 
modes encountered when we quantize a neutral field. Indeed for the charged field the two 
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linearly independent solutions of the radial equation ( f4.2| ) have to be taken into account 
while for the uncharged field only one has to be considered, the other blowing up exponen- 
tially at infinity. Here again the analysis of the classical trajectories helps to understand 
what happens. The modes that have to be rejected correspond to particles pushed to 
spatial infinity, i.e. quantum mechanically they describe fluctuations engendered by the 
Schwinger mechanism and disappear when E t— > oo. 

This indicates for the accelerated detector as well as for the charged black hole, that the 
Schwinger process first quickly switches off the external field, by emission of preferentially 
charged particles that will neutralize the plates of the condenser producing the electric 
field, or the black hole. Only then will the system reach a thermal (quasi)-equilibrium 
state at the Unruh/ Hawking temperature. 

• The ratio of population of a heavy two level charged detector accelerated and in interaction 
with the charged scalar field are in accordance with thermodynamics, as far as the charge 
difference between the two levels of the detector is small compared to their difference in 
masses. Moreover the amplitude of transition of such a detector can be obtained as the 
large mass limit of a model built on three interacting fields, a model where the levels of the 
detector are quantized. This shows that it is only in the limit where the recoil effects and 
pair creation a la Schwinger of the quanta describing the detector levels become negligible 
that the thermodynamical limit makes sense and is recovered. It is also noteworthy that, 
as expected, in the limit of electric field going to zero, the probabilities of transition by 
absorption or emission ( |7.21| , [7.23| ) of the fluctuating modes whose distributions become 
null ( |6.32|) also vanish. 
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A Asymptotic behavior of Whittaker's functions 



The Rindler's modes are built from the solutions of eq.( |4.2|) , which are given by Whittaker's 
functions. For completeness, we recall in this appendix the few properties of these functions 
used in the main text. 

Defining W + = exp(a^)JF as a function of the imaginary variable z = — ie-§^ exp(2a£), it is 
immediate to verify that W + [z] has to satisfy the Whittaker's equation: 



dz 2 



+ 



1 K 1 - 4/i 2 

I h — 

4 z Az 2 



W 







with 



% l lo m \ uj 

k = — , u = i — 

2\a qEJ ' p 2a 

The Whittaker's functions are related to the confluent hypergeometric function by 



W + [z] = e- z/2 z 1/2+ »F{K^,z] 

where F[k,/i, z] is a solution of the Kummer's equation: 

id 2 d 1 1 

I Z dz^ + ( 1 + 2 V- z )-j~ ~ ^~ K ) \ F K V, A 
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(A.2) 



(A.3) 



(A.4) 



In momentum space (p = —d/dz), this equation becomes first order. Solving it by an elemen- 
tary quadrature and returning to the z representation, we obtain an integral representation of 
the Whittaker's function: 
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(A.5) 



The second of these equations results from the change of variable q = pz and the Jordan's 
lemma allowing to integrate on the real axis instead of the imaginary one. When z goes to 
infinity, eq. (|A.5 ) gives immediately the asymptotic behavior of the function: 



W + [z] ~ e~ z/2 z K 



(A.6) 



while near z = 0, the Whittaker's function behaves like 

r[2/i] zl/2 ^ + r[-2//] 



W + \z] ~ 

L 1 z->0 



r[± + fi-K 



r[| -fl-K] 



(A.7) 



If, instead of the change of function ( |A.3|) , we use the one with the opposite sign of /i, setting: 

W+[z] = e~ z/2 z 1/2 -^F[K, -n, z\ (A.8) 
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we obtain a different integral representation, but of the same function (see eq.( |A.7| )): 



r[| — ji — k] Jo 



r°° i i 

/ e - pz P -2-»- K (i+p)-2-» +K dp = w+y 

Jo 



(A.9) 



Indeed, both representations are solutions of a second order differential equation and have the 
same asymptotic behaviors. In the mathematical literature, this function is usually denoted 
by: "W K ,(i(z), an d satisfies the relation : 

W + [z] = W K>+fl (z) = W K ,_ M (z) . (A. 10) 

A second, independent, solution of eq. QA.l ), is given by the complex conjugate of the first one 

W~[z] = \W + [z}]* . (A.ll) 
So, we obtain complete sets of unnormalized modes, solutions of the wave equation (|4.1| ): 
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(A.12) 



From the asymptotic expansions of the function W*, we read immediately the various coefficient 
that define the charges carried by these modes: 
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(A.13) 
(A.14) 

(A.15) 
(A.16) 
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In the main text, we also make use of linear combinations of the previous Whittaker's func- 
tion whose behaviors near £ = — oo, [z = 0] are particularly simple. They are the functions 
M± Kjfl (±z) and its complex conjugate M TK _^(^fz): 



\aj a 



e 2 i E 



r 


'1 _ 






to 





2 + * 2g£;)j 



w 



obeying the relations 



M K Jz) 



M K ,-»(z) 



■ tt I TV omega 

e~^e +e ^^M_ K J-z) 



K-fl 



from which we have defined the (unnormalized) modes: 
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whose charge content is given by the coefficients: 
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(A.21) 



(A.22) 

(A.23) 
(A.24) 



(A.25) 
(A.26) 



(A.27) 



sgnH \D + (Mt)\ 2 = e (\C + {Mf)\ 2 - \C.{Mt) 

Let us remark that W+(r,f) = [W+(-T,f)]* and A4 + (r,0 = [-M+(-r,0]*, i.e. in 
classes of modes are related by a r inversion followed by a complex conjugation. 



(A.28) 

(A.29) 
and out 



B Explicit forms of some Bogoljubov coefficients 

For sake of completeness we recall here the basic relations between Fock basis defined trough 
Bogoljubov transformations and give here the explicit expressions of the various Bogoljubov 
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coefficients that we have used in the main text. Suppose a quantum field operator satisfying 
a second order field equation, is defined in terms of two sets of linearly independent solutions 
(pi, 4>l* and $^*, labeled by the indices k and K, of the corresponding classical equation. 
We have 



* = n^i+m* (b.i) 

k 

= Y. A k^k + B ] k ^ k * , (B.2) 

K 



and 



In* = E^i^+^r Wr = ^t^r-pt^c [ } 

If, as usual, we suppose these basis orthonormal, unitarity implies : 

W = Ei4*^-7A*7*. 1 = Eje/* e / (B.4) 



Combining these relations we obtain the links between the various creation and annihilation 
operators : 

(S = 5^ a i + oW and similarlly Mf = ^Tt)~fi«] ' 

(B.5) 

For each set of operator is associate a "vacuum " state \ fl > and \uj > such that for all quantum 
number k and K : 

A K \Vl >= B K \Q >= and a k \co >= b k \co >= . (B.6) 
A standard computation |30], [8| gives the link between these vacuum states : 

\n >= a/VE,,™^} \ u> , \uj >= a/VE^ kl 4X} |n > (b.7) 

where the matrices m kl and are given by : 

ro « = 53(a- 1 )}A and M M = EriK 1 )^ • (B.8) 

The normalization coefficients are of particular physical interest; they give the projection of 
one vacuum state on the other (a fact that we anticipate in the notationJV*) . In the main 
text we only encounter "diagonal" Bogoljubov transformations, i.e. such that the two sets of 
indices {k} and {K} are identical and the various matrices a k Kl ... e|- oc 8\ i.e. diagonal. 
In such cases the matrices m kl and M KL are also diagonal and the normalization coefficient 
are particularly easy to evaluate. By factorizing the vacuum state according to the quantum 
number k: \uo >= Y\ k \uo k >, we obtain : 

00 1 , 

1 " ■ ,k|2 
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a 



k I 



(B.9) 

the phase of M remaining arbitrary, the vacuum state being actually ray in the Hilbert space. 
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Bogoljubov transformations between Rindler in and out modes 

On quadrant R , the coefficients connecting Rindler in and out-modes, coefficients occurring 
in eq.( |4.26|) , are given by: 



> 0) = omec/a < 0) 



, _™ r[<*] r[|-z(^-fi)]^(^) 



" r[-ia] r[i + z^] AT(vr^) ' 
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On quadrant P , see eq.( |4.54j ) they are given by: 
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Bogoljubov transformation between Unruh and Rindler out modes 

In the R and L sectors, Unruh out modes are related to Rindler modes as: 
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TfflkM = e$>) = a$;M=/^>). (B.13) 

In the same way, we get the expression of the Unruh out modes in terms of Rindler P and F 
modes 

wl out = flM,(a))^ + M«Su,W (B.15) 
= ^»V™1 + {e$»V-j> + 78?>)Z<#} + 

{t8S»K? + ^(^W} (B.16) 

with coefficients: 
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(w > 0) = a8&, (a, < 0) = 7 St* (w > 0) = e^ p * (w < 0) = (w < 0) 
/?X > 0) = «1 (w < 0) = C*; (w < 0) = 7 X (w > 0) = /4" v * p (a; < 0). 

(B.18) 

C Schwinger representation for wave function products 



In subsection |5.1| we made use of several integral representations of products of modes. Here- 
after, we briefly indicate how they are obtained. 



•Inertial observer Using the relations (9.240) and (6.669.3) of ref. [13|, we may express the 
product of modes: ip^. out (x)(p^ out (x) as: 



■ 2 w m+i^m+i^) 

r» ^L_ e ^(^^)(-^-sinho (coth i rl f|_ (c i; 

Jo (sinh^) _ 2 2 

Introducing the new variable s defined by sinh2gi?s = (sinh£) _1 in this integral, we get: 
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(C.2) 

•Accelerated observer In the same way, the product of wave functions U™ t (x)V™*(x) can be 
expressed as the integral : 

UTm - - i f,^ T j_ r( , +i (^g); +ia) 
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r jtM-frt&^i^j??) (c.3) 



thanks to the eq. (6.669.6) of [13 . 
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Figures captions 



Fig. 1 The four Rindler patches and their coordinates. 

Fig. 2 Typical trajectories of charged particles on Rindler patches: plain curves correspond to 
particles trajectories, the dashed one to an antiparticle trjectory. The numbers labelling 
the different trajectories refer to the discussion in the main text. 

Fig. 3 Schematic representation, in the various Rindler quadrant, of the charged carried by 
the various modes defined in section |3|, and of typical wavepackets built out of them 
with the appropriate mean value of uj. On the R and L quadrants, we consider packets 
centered around a positive mean lo\ on the P and F quadrants, the packets are built out 
of modes whose u are essentially negative. 

Fig. 4 In and out Unruh modes represented as superposition of Rindler modes and their 
charge contents 

Fig. 5 Ratio of the total probabilities of excitation and deexcitation of a charged detector of 
mass gap AM /a = 0.1, normalized with the Boltzman factor exp — [27ru;*/a], as function 
of the mass m/a and "charge" qE / a 2 of the exchanged quantum. 
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